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Abstract. We prove the existence of the double scaling limit in the unitary matrix model 
with quartic interaction, and we show that the correlation functions in the double scaling limit are 
expressed in terms of the integrable kernel determined by the psi-function for the Hastings-McLeod 
solution to the Painleve II equation. The proof is based on the Riemann-Hilbert approach. 
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1. Introduction 



In this paper we will concern with the double scaling limit in the unitary random matrix model 
with quartic interaction. The unitary random matrix model, or the unitary ensemble of random 
matrices, is defined by the probability distribution 

^NidM) = Z-^ cxp{-NTrV{M))dM, Zn = [ exp {-NTrV{M)) dM, (1.1) 

on the space Hn of Hermitian N x N matrices M = (-^ij)i<j j<jv> where in general V{M) is a 
polynomial of even degree with a positive leading coefficient, or even more generally, a real analytic 
function with some conditions at infinity. The basic case for the double scaling limit is the quartic 
matrix model when 

V{M) = ^M^ + |m^, g>0, (1.2) 

and we will consider this case only. By a change of variable one can reduce the general case to 
the one with g = 1 but we prefer to keep g because it is useful in some questions. The double 
scaling limit describes the asymptotics of correlation functions between eigenvalues in the limit 
when simultaneously iV — > oo and t approaches the critical value 

tc = -2^5, 

with an appropriate relation between N and t — t^. (see below). The critical value tc is a bifurcation 
point: for t > tc the support of the limiting distribution of eigenvalues consists of one interval, 
while for t < tc it consists of two intervals. This can be described as follows. 

Let du{x) be the limiting eigenvalue distribution (for the proof of the existence of the limit- 
ing eigenvalue distribution and formulae for it see [BPS], [DKM] and [DKMVZl]; see also earlier 
physical works [BIPZ], [BIZ], and others). The distribution diy{x) is a unique solution to a vari- 
ational problem and for a general analytic function V{M) satisfying certain conditions at infinity 
(see [DKM] and [KM] for more detail) the distribution di'{x) is supported by a finite number of 
segments [oi, 6i], . . . , [a^, bg], it is absolutely continuous with respect to dx, and on J = ^'j=i[0'j,bj] 
its density function is of the form 

P{x) = —hix)R]/'^{x), R{x) = Y[{x - aj)(x - bj), 

1/2 

where h{x) is a polynomial and (x) means the value on the upper cut of the principal sheet of 
the function i2^/^(z) with cuts on J. The variational problem for du{x) implies that the polynomial 
h{z) satisfies the equation 

^ ^ - (1.3) 

where 

^ r p{x)dx ^ ^ ^ 2)^ ^ ^ ^ 

J. J z-x 

(see e.g. [DKMVZl]). Equation (1.3) enables one (sec e.g. [BPS]) to determine uniquely the limiting 
eigenvalue distribution for quartic polynomial (1.2). For t > tc the support of the distribution 
consists of one interval [— a,o] and on [—a, a], 

p(x) = -Cbo + b2X^)\Jo?- - a;2, \x\ < a; t> tc- (1.4) 

TT 
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Here 



a = 



-2t + (4*2 + 48^) 



1/2 N 1/2 



35 



g t+((tV4)+3ff) 
b2-^, bo- 3 



1/2 



For t = tc, (1.4) reduces to 



p{x) 



27r 



(1.5) 



where a = l/g'^/^. Observe that f»(0) = in this case. For t < the support consists of two 
intervals, [—a, —h\ and a], and on these intervals, 



1 



p{x) = -bQ\x\^/{a? - a;2)(a;2 - h^)^ < 6 < |a;| < a; t<tc. 

TT 



Here 



-2^5 



1/2 



6o = 2- 



5 / V 5 

Fig.l depicts the density function p{x) for 5 = 1 and t = —1, —2, —3. In this case tc = —2. 



(1.6) 



(1.7) 




Fig 1: The density function p{x) for 5 = 1 and t = —1, —2, —3. 



Correlation Functions and Orthogonal Polynomials. Our basic object of interest is correlations 
between eigenvalues. The m-point correlation function KN^ni^i, ■ ■ ■, Zm), m = 1, 2, . . . , is a dis- 
tribution over the space X'(R™) of C°°-functions ip{zi, . . . , Zm) on with compact support such 
that for the product functions (p{z-\_, . . . , Zm) = fiizi) ■ ■ ■ <Pm{zm), 

KNmiM^l)...<PmiZm))= [ [x^, (M) . . . X,,^ (M)] /Xjv (c?M) (1.8) 



where 



N 



Xv.(^) = IZ^(Ai), Mej = \jej. 
i=i 

By linearity K^rni'^) extended then from the set of product functions ipi{zi) . . . fmizm) to the 
whole space V{W^). 

The m-point correlation function K]yyn{zi, . . . ,Zn) turns out to be a regular function on the 
set {zi / Zj} and on this set it can be expressed in terms of orthogonal polynomials. Namely, let 



P„(z) = 2;" + ..., n = 0,1,2,..., 



(1.9) 
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be monic orthogonal polynomials on the line with respect to the weight e ^^(^) , so that 

Pn{z)Pm{z) e-^^(^)dz = hnSmn. (1.10) 



/ 



We normalize Pn{z) by the condition that the leading coefficient of Pn{z) is equal to 1. The 
polynomials Pn{z) satisfy the recursive equation [Sze] 



zPniz) = Pn+l{z) + RnPn-l{z), (1.11) 

where 

/?„ = 

hn-l 

We associate with P„(z) the functions 



= > 0. (1.12) 



^„(z) = ^P„(^)e-^^(^)/2, (1.13) 

\/hr, 



'•n 



which form an orthonormal basis in L 

/oo 
tpniz)tpm{z) dz = Snm- (1-14) 
-oo 

Prom (1.11), 



Z^Z) = i?y+Vn+l(-^) + Rl/^^n-l{z), (1.15) 

The m-point correlation function K^^^zi, . . . , 2;„) is expressed in terms of the functions ipn{z) as 

KNm{zi,...,Zm) =det{QN{Zi,Zj)).j^-^^ Zi^Zj, (1.16) 

where 

N-l 

Qn{z,w) = ^ '>pk{z)'ilJk{w) (1.17) 

fe=0 

([Dys]; see also [Meh], [TWl]). The Christoffel-Darboux identity reduces Qn{z,w) to 



n pi/2 ^n{z)^n-i{w) - i;N-i{z)ipN{w) 
LJn{z,w) = Kj^ ■ l-L-J-oj 

z — w 



We are interested in the asymptotics of the correlation functions K^mizi, . . . , Zn) in the limit 
when iV — >■ oo. Formula (1.16) reduces the problem to the asymptotics of the kernel Qn{z,w) 
and (1.18) further to the asymptotics of the functions tpN{z) and ipN-i{z). In noncritical case 
t < tc this problem was solved in [BI]. Noncritical asymptotics for the general V{M) were found in 
[DKMVZl]. 

Double scaling limit for correlation functions. We will consider the correlation functions 
KNm{zi, ■ ■ ■ , Zm',t) in the situation when the parameter t approaches the critical value tc and 
Zi, . . . ,Zm are near 0. The problem is how to scale eigenvalues and t — tc to get a nontrivial limit 
for the correlation functions. More precisely, we want to find numbers ^ and ij, critical exponents, 
such that the following "double scaling" limit exists: 



/Z^ Zm . 



lim 



;te + yiV-«) = Km{zi,...,zm;y), (1.19) 



and it is a nonconstant function of the parameter y, — oo < y < oo. We also want to find 
the hmiting correlation functions Kmizi, . . . , Zjn',y)- We will derive the double scaling limit for 
correlation functions from corresponding asymptotics for orthogonal polynomials. The relation 
between the double scaling limit for correlation functions and the one for orthogonal polynomials 
can be described as follows. 

Double scaling limit for orthogonal polynomials. Let us fix the parameters t < and g > 
in quartic polynomial (1.2). Since t < 0, V{z) is a double- well function. Consider orthogonal 
polynomials Pn{z) with respect to the weight e~^^^^\ The bifurcation of p{x) in the quartic 
matrix model is closely related to the bifurcation of the distribution of zeros of the orthogonal 
polynomials Pn{z)- This relation is motivated by the Heine formula for orthogonal polynomials as 
follows. Let xi,. . . , x„ be zeros of Pn{z) which are all real. Consider the corresponding probability 
measure 



dfiri{x) = n ^ 6{x — Xj) dx. 



We can write the weight as 



g-w(x) ^ g-nv.(x)^ y^(^) _ A-iy(x), A = 
By the Heine formula for orthogonal polynomials (see e.g. [BIZ]), 

P„(z) = (det(z - M))y^ ^ Z-\Vx) J det{z - M)e-"^^^(^)dM, 
or equivalently, in the ensemble of eigenvalues, 

Pn{z) = (fl{z-XA . 

Due to the variational principle, we expect that if n is large then the probability measure 

n 

du{x; {Xj}) = 6{x — Xj)dX 

i=i 

is close to the equilibrium measure p{x; Vx)dx for typical {Xj}. Therefore, 

n"^log( JJ(2;- A,) \ log{z - x)p{x;Vx)dx, 

hence by Heine, 

/\og{z — x)d^Xn{x) = log Pn{z) K, I \og{z — x)p{x\V\)dx, 2:GC\M. 
-oo J J\ 

That is, as n — >■ oo the logarithmic potential of diJ,n{x) converges to the logarithmic potential of 
p{x] V\)dx. Thus, we can expect that as n, AT — > oo in such a way that ^ A, the measure 
diJLn{x) converges to d/Xoo(x;A) = p{x\Vx)dx. This convergence has been rigorously established 
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for the general V{M) in [DKMVZl] (see also [BPS], [BI] and works on the theory of orthogonal 
polynomials [LS], [Lub], [Nev], [ST], and references therein). 

For A small, the zeros of Pniz) are located near the minima of the function V{z) and the 
support of d^oo{x\ A) consists of two intervals. For A big, the support of dfioo{x; A) consists of one 
interval, and there exists a critical value, 



(1.20) 



which separates the two- interval and one-interval regimes (see [BI]). The problem of double scaling 
limit for orthogonal polynomials is to prove the existence of a nontrivial scaling limit for the tp- 
function. 



lim 

n,JV->oo: (n/iV)=Ae+yJV-« 



with some critical exponents ^ and r], where Cn 7^ are some normalizing constants. Our main 
goal in this paper is to derive uniform asymptotics for the functions tpni^) on the whole complex 
plane as n, A?^ ^ cxd in the double scaling limit regime, (n/N) = Ac + yN~^. When applied to the 
quartic matrix model, these asymptotics will give us the double scaling limit for the correlation 
functions at the origin as well as scaling limits in the bulk of the spectrum and at the edges. There 
are three basic ingredients in our approach: the string (Freud) equation, the Lax pair for the string 
equation, and the Riemann-Hilbert problem. 

String Equation. Let Q = iQmn)m,n=o,i,... be the matrix of the operator of multiplication by 
z in the basis ipniz)- By (1.15), 



Q 



( 



\ 



o_ 














Let P be the matrix of the operator d/dz. Then P* 
Observe that 



...X 

-P, where P* is the transposed matrix. 



pi p-NV/2 

= -{NV'/2)^^ + " ^ = -(W/2)V'„ + 



'hr. 



n 



n-l 



+ 



SO that 



and 



hence 



[P + NV'{Q)/2\ 



n 



n,n— 1 



'Rn 



(1.21) 



0=[P + W(Q)/2]„_,^„ 
Combining (1.21) with (1.22) we obtain that 

n 



[P + NV'{Q)/2]^^^^, = Q, 

-P + NV'{Q)/2l^^_, 



njr: 



= [v'mn.. 



n—lj 



(1.22) 



(1.23) 



the discrete string equation (see [BK], [DS], and [GM]). When V is as in (1.2), this reduces to 
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n 



— Rn{t + gRn-1 + gRn + gRn+i) ■ 

(cf. [Fre], [BIPZ], [IZ], and [BIZ]). Note that from (1.12) and (1.24) it follows that 

— < tRn + gRn, 



hence 



< i?„ < 



25 



Lax Pair. Indroduce the vector valued function 



V'n-l(^) ) ' 



It satisfies the following 2x2 matrix diflFerential equation: 
where 



-Ci + ^+gzRn) 
-VR^{t + gz^ + gRa-i + gRn) 



+ gz^ + gRn + gRn+i) 
'-i + ^+gzRn 



(1.24) 



(1.25) 



(1.26) 



(1.27) 



(1.28) 



([FIK2]; see also [BI]). For the sake of completeness, let us remark, that for a general even polyno- 
mial 



- 2, 



the matrix An{z) has the form 

Aniz) = 



' Rj 



-an{z) 
■ an-i{z) + an{z) 



I Rr 



Uniz) + an+liz) \ 



an{z) 



where 



aniz) = 



V'iz) 



J=2 1 = 1 



(see [FIK4]). This reduces to (1.28) when V is as in (1.2). 
In addition, we have the recurrence equation 



where 



= Un{z)^n{z), 



Un{z) 



pl/2 „l/2 
1 







(1.29) 
(1.30) 



The compatibility condition of equations (1.27) and (1.29) is 

C/;(z) = NAn+iiz)Un{z) - NUniz)Aniz). (1.31) 

Restricting this matrix equation to the element 11 one can rederive string equation (1.24) (see 
[FIK]). And vice versa, (1.24) implies (1.31). Thus, system (1.27,29) is the Lax pair for equation 
(1.24). 

Remark. The Lax pair (1.27,29) can be alternatively written (see [FIK]) as the following linear 
differential-difference system, 

dz , V = ('01,^2, ■.■,V'n,.-.). 

Qip = zip, 

Accordingly, string equation (1.42) takes the "quantum" (cf. [Moo], [Nov]) commutator form, 

[Q,vLm = ^ 

(see [FIKl]). Here, as usual, the low triangular part, M_ , of a matrix M is defined by the equations, 

[M_]nm = [M]nm if n > m and [M_]„„ = if n < m. 

To describe the Riemann-Hilbert problem consider the adjoint functions 



4>niz) = e 



r e-'"-'->/V„M ,,c\M. (1.32) 

7_oo U-Z 



They share the following properties: 

(1) As Z — GO, 

oo 

72fc 

(2) The function (j)n{z) has limits 4'n±{z) as z approaches the real axis from above and below, and 
the limits are related as 

(3) The vector function 

(t)n{z) 



^n{z) 



(t>n-l{z) 



solves Lax pair equations (1.27,29). 

To prove (1) expand into geometric series and use orthogonality. Property (2) follows from 
the jump condition for Cauchy integral. Property (3) will follow from the Riemann-Hilbert problem 
(see Proposition 1.1 below). Define the matrix- valued function *„(z) on the complex plane as 

^n{z)=( f-^'), /"^'\V (1.33) 

It solves the Lax pair equations, 

U'^{z)=NAr,{z)^r.{z), 

I *„+i(z) = C/„(z)*„(z), ^ ■ ^ 



and the Riemann-Hilbert problem for orthogonal polynomials. 

Riemann-Hilbert Problem. The Riemann-Hilbert problem for orthogonal polynomials is for- 
mulated as follows ([FIK]; see also [BI], and [DKMVZl]). One has to find a 2 x 2 matrix- valued 
function '^n{z) on the complex plane which is analytic outside of the real line, and which has 
continuous limits from above and below of the real line, 

'^n±{z)= lim ^^n{u), 

u—>z, ±lmu>0 

SO that 

(i) "ifniz) satisfies the jump condition on the real line, 

^n+{z) = ^n-{z)S, Rez = 0, (1.35) 

where 

s={l -f); (1.36) 

and 

(ii) As z — oo, the function ^^(z) has the following uniform asymptotic expansion: 

where T/., A; = 0, 1, 2, ... , are some constant 2x2 matrices, with 

Xn is a constant, and as is the Pauli matrix, 

'3= (J _°i). (1.39) 

To solve the Riemann-Hilbert problem we have to find the function ^ni^) real numbers 
Xn and Rn > such that (i) and (ii) hold. 

Proposition 1.1. There exists a unique solution to the Riemann-Hilbert problem (1.35-38), 
and it is given by (1.33). The number Rn in this solution coincides with the recursive coefficient 
in (1-11) and 

K^"^. (1.40) 

where hn is defined in (1.10). The functions ^'„(z), n = 0, 1, 2, . . . , satisfy the Lax pair equations 
(1.34). 

For the proof of Proposition 1.1 sec [FIK], [BI], and [DKMVZ2]. From now on, we shall "forget" 
(cf. [BI]) about explicit equation (1.33) for the solution ^^(z) of the Riemann-Hilbert problem 
(1.35-38). It is the Riemann-Hilbert problem itself which becomes now the principal characteristic 

of the function ^'„(z). In other words, we shift the focus from equation (1.33), which represents 
the function ^n{z) in terms of the orthogonal polynomials Pn{z), to the equation 



V'n(^) = [*n(^)]ll, 
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which represents the orthogonal polynomials Pn{z) in terms of the function ^„(z) (which in turn 
is uniquely determined as a solution of the Riemann-Hilbert problem (1.35-38)). It also worth to 
emphasize that (see [BI]) in the setting of the Riemann-Hilbert problem the quantities Rn and 
A„ are not the given data. They are evaluated via the solution *„(2;), which is determined by 
conditions (1.35-38) uniquely without any prior specification of Rn and A„. 
Since V{—z) = V{z), we obtain that 

^ni-z) = (-l)"Vn(^)- (1-41) 

This leads to the following equation on '^ri{z): 

^ni-z) = {-iras^niz)<73. (1.42) 

We turn now to formulation of the main result concerning the semiclassical asymptotics of the 
functions "ipn{z) in the double scaling limit. 

Formulation of the main result. We will consider n such that n/N is close to the critical value 
Ac = t^ /{Ag). We introduce the variable y as 

V = ^'^'''{^-^^). <*=(£)'"• (!•«) 

We will assume that \y\ is bounded, that is we fix an arbitrary large number Tq and we will consider 
such n that 

\y\ < To. (1.44) 

Our results will concern asymptotics of the recurrence coefficients Rn and the functions ipniz). For 
Rn we will prove the following asymptotics: 

Rn=^^+ N-'/^c^i-ir+Mv) + N-^^^C2v{y) + O(Ar-i), 



2\t\ \ 1 / 1 \ 1/3 



(1.45) 



where u{y) is the Hastings-McLeod solution to the Painleve II equation 

u"{y) = yu{y) + 2u^{y), (1.46) 
which is characterized by the conditions at infinity, 

lim =1, limu(y) = 0, (1.47) 

and 

v{y) = y + 2u\y). (1.48) 

Remark. Asymptotics (1.45) has been suggested in physical papers by Douglas, Seiberg, 
Shenker [DSS], Crnkovic, Moore [CM], and Periwal, Shevitz [PeS]. The existence and uniqueness of 
the solution to (1.46,47) was first established in [HM] (see also the later works [Kapl,2] and [DZ2]). 
It is also worth noticing that the asymptotics of u{y) at -|-oo can be specified as 

u(y)= Ai(y)(l + 0(H-i)), (1.47') 
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where Ai (y) is the Airy function. Moreover, the asymptotic condition (1.47') characterizes solution 
u{x) uniquely, so that the first equation in (1.47) and equation (1.47') constitute an example of 
the so-called connection formulae for Painleve equations (see [IN], [FI], and [I] for more on this 
matter). We should also point out that the analytic and asymptotic properties of the Hastings- 
McLeod solution for complex values of y can be extracted from the results of [Kapl,2], [Novok], 
[Kit2], and [IK]. 

We will prove asymptotics (1.45) simultaneously with semiclassical asymptotics for the func- 
tion ijjn{z) on the whole complex plane. To that end we will substitute asymptotics (1.45) into 
the coefficients of differential equation (1.27) and solve the resulting equation in the semiclassical 
approximation. Realization of this idea is not trivial and in comparison to the standard semiclas- 
sical analysis, it involves a new type of special functions generated by the monodromy problem 
associated with the second Painleve equation. The appearance of the latter is closely related to the 
fact that ansatz (1.45) leads to the coalescence of four turning points in system (1.27) (cf. [Kitl]). 
To see this let us analyze (1-27). 

System (1.27) can be reduced to one equation of the second order. Namely, if we denote by 
(^ij J hj = 1) 2, the matrix elements of the matrix An{z), solve V'n-i in terms of from the first 
equation in (1.27), 

V'n-i = (A^~ Vn - ailV'n) , 



and substitute 



Ol2 



(1.49) 



then we obtain from the second equation in (1.27) the Schrodinger equation on 77, 



(1.50) 



where 



[/ = -det^„ + AT-i 



(an)' - an 



ai2 



- N- 



(ai2)" 3((ai2)') 



l\2-\ 



2a 



12 



4(ai2)2 



(1.51) 



The zeroth order term in the potential U is — det An so let us consider properties of — det 
Prom (1.28) we obtain that — det An{z) is a polynomial of the sixth degree, 

2 4 

- det An{z) = ^{z^ - zi) + a„z2 + /3„, 



where 
and 



zo 



an = ^ - gRn{t + gRn-1 + gRn + gRn+i), 

Pn = —RnOn-lOn, On = t + gRn + gRn+1- 

Due to string equation (1.24), a„ simplifies to 



(1.52) 
(1.53) 
(1.54) 



gn 



an = ^-'-TT = -g[^- \c) , 



(1.55) 



or, according to scaling (1.43), 



an = -gcoyN ^/s. 
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The substitution of asymptotics (1.45) into /?„ gives 
where 

'^l/3|i|l/3' 



C3 



25/3 



-j[v^{y)-Aw\y))], C4 = (-ir(^,f^ 



w{y) 



and ■«;(?/) = Thus, 



^2 ^4 



— det^„(2;) 



As iV — oo, — det A„(2;) approaches the polynomial 

4 



Oooi^) = — i-(2 -^o)> 



(1.56) 
(1.57) 

(1.58) 

(1.59) 



which has two simple roots at ±zo and a quadruple root at the origin. Therefore, equation (1.50) 
has two simple turning points approaching ztz^ and four turning points coalescing at the origin. 

In accord with this analysis, we divide the complex plane into several regions. The form of 
the semiclassical asymptotics will be different in different regions. Let di and d2 be arbitrary fixed 
(i.e. independent of N) numbers such that 



Zn 

0<d2<di< ^. 

4 

Introduce the rectangular region CI as 

Q = {z : |Re z\ < zq + di, |Im z\ < ^2}. 



(1.60) 



(1.61) 



Im z 




Fig 2: The region J7. 



Observe that the segment [— zo,-2o] lies in il. Introduce furthemore the rectangular subregions 
and of the region Q as 

n^ = {z: \Rez\<di, |Imz| < 0^2}, = {z : \Rez- zo\<di, |Imz|<d2}, (1.62) 

so that Q'^ is a rectangular region centered at and at zq. Finally, let J7i be the rectangular 
region between $7° and Cl^, that is 

fti = {z : di < Re z < Zq — di, \lmz\ < ^2} (1.63) 
12 



(sec Fig. 2). We will prove the following semiclassical asymptotics for the vector function ^'^(z): 

1. In = C \ $7: WKB asymptotics of exponential type. 

2. In Oi: WKB asymptotics of cosine type. 

3. In Q^: turning point (TP) asymptotics (in terms of the Airy function). 

4. In ri°: critical point (CP) asymptotics of Painleve II type. 

In the regions —il^ and — r^i the semiclassical asymptotics will follow by symmetry equation (1.42). 
As a matter of fact, both the TP and CP asymptotics will be extended to the region f^i and we 
will use this exte nsion t o prove the WKB asymptotics of the cosine type in fli. 

Region = C\Cl, WKB asymptotics. In $7^^ we will prove the following asymptotics: 



l + O 



N{l + \z\ 



WKB 



where 



-iVr(2) 

-ivr(z) 



with the quantities on the right defined as follows. We define 

= 1^ + iV-V3ci(-l)"+i„(y) + N-^/^C2v{y) 
[cf. (1.45)], where y is defined as in (1.43), and we put 



(1.64) 
(1.65) 

(1.66) 



tz gz 



+ gzR\ 



^^{t + gz^ + gRl + gRl^^ 



(1.67) 



2 2 

\-^^{t + gz'' + gRl_^+gRl) 
[cf. (1.28)]. We denote the matrix elements of A^n{z) by aPij{z). We define the function //'^(z) as 



where U'^{z) is a suitable approximation of the function U{z) in (1.51), namely. 



(1.68) 



and we put 



du. 



(1.69) 



where Zq is the root of U^{z) that approaches zq as A" — oo. The contour of integration in (1.69) 
is taken as follows: first from Zq to zq + di and then from Zq + di to z around the region Q, in the 
counterclockwise direction. The matrix T'^{z) is defined as 



li-{z) 



1/2 



\ 



ali(z) 



(1.70) 
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Observe that detr^(z) = 1. The branches for (a?2(^))^/^, (^7°(^))^/^ and (^=(z))i/2 are fixed by 
the condition that they are positive for large positive z. We will check below that the function 
^wkb(2;) is analytic in O'^. The meaning of formula (1.64) and similar formulae to follow is that 
there exists a 2 x 2 matrix valued function ej^{z) such that 



and 



= (l+£iv(z))^'WKB(2), Z^9F, 



sup 1(1 + |z|)£jv(-z)| < cA^"\ c>0. 



For concrete calculations the function \j'^{z) can be simplified as follows. The function 

[/!(.) = a?/(.)-a?,(.)^'^") 



is analytic in OF and as AT — oo, 



23 = (-l)"2i/3^i/3|i|i/3^ ~ 



22/3g2/3 [y(^y)+4{-irw{y)] 
4|t|i/3 



(1.71) 



hence 



,2 ^4 



t/°(z) = ^(.^-.o^)- 



n 



/3iv = (C3 - cs)N-^/^ + (C4 - C4)Ar 



-5/3 



and the error term 0{N~'^) in U^(^z) can be neglected in concrete calculations, reducing U^(z) to 
a polynomial of the sixth degree. 

For the normalizing constant /i„ we will prove the asymptotics 



/in = e -^^o" (l + 0(iV-i)). 



(1.72) 



The integral in the exponent is regularized at infinity as in formula (5.3) below. 

Region Qi, WKB asymptotics of cosine type. In the region J7i we define ^vvfKB(^) as 



*wkb(^) = ^«)-'/'Ti(^) 



( cos(iVei(z) + ^) ^ 



sm 



z G rii, 



(1.73) 



where 



and 



f,^{z) = {-U\z)f'\ ^^{z)= !\^{u)du- zeni, 



Ti{z) 




(1.74) 



(1.75) 
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so that detTi(z) = 1. The functions —U^{z) and ai2iz) are positive on the interval di < z < ZQ — di 
and the branches for the square roots in (1.74,75) are chosen to be positive on this intervaL We 
will prove that 

*„(z) = (l + 0(iV-i))*WKB(^), zeni. (1.76) 

Region il^, TP asymptotics. In the region we define the turning point vector valued function 
^TP (z) as 

where 

^i^) = (I II VU^) duj . (1.78) 

and 

/ ^ \ \ 1/2 / 1 \ 

W{z)= ^44) ^ , (1.79) 



w'{z) 



ai2(^) «i2(^) 



so that det W{z) = 1. The branches for fractional powers in (1.78,79) are chosen to be positive for 
z > Zq . Observe that the function w{z) is analytic and has no critical points in fl^. We will prove 
that 

^niz) = {l + 0{N-'))^Tp{z), zGn\ (1.80) 

Region 0°, CP asymptotics. The normal form for system (1.27) at the critical point 2; = is 
the system 

*'(z) = A{z)^{z), (1.81) 



where 



4 . ^ _ / (-l)Mn(y)z + (-l)-2«;(y) + v{y) , 

- 1 -4^2 + {-l)^2w{y) - v{y)) -{-l)''4u{y)z > ' ^ ' 



and 

z;(y) = y + 2t/2(y), wiy) = u'iy), (1.83) 

where n(y) is the Hastings-McLeod solution of the Painleve II equation (1.46) (see Section 3 below). 
We will consider a special solution to (1.81), 

*W = (|^[Ji). d-s") 

which is characterized by the following properties: 
1) ^{z) is real, i.e.. 



$(z) = $(z). (1.85) 

2) It satisfies the parity equation 

i{-z) = (-l)"(73$(z). (1.86) 

3) On the real axis the functions {z) have the asymptotics 



$i(z) = cos ( — + - — j + O (z-i) , 

$2(z) = -sm( — J +0(2 ^), z^ioo. 



(1.87) 
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The existence of the solution ^(z) is a nontrivial fact of modern theory of Painleve equations (see 
e.g. [IN], [DZ2]). It should be also noticed that, in addition to properties (1.85-87), the function 
$(z) is an entire function on the complex z plane, and its asymptotic behavior is known in the 
whole neighborhood of 2; = oo (see Proposition 3.2 below). 

Remark. As a function of the parameter y, the vector ^{z) is a meromorphic function which 
satisfies the linear differential equation (cf. (1.81)), 

^ = B{z)m, (1-81') 

where 

Moreover, the asymptotic distributions of the complex poles (which coincide with the poles of the 
Painleve function u{y)) and the large |y| asymptotics of ^{z) can be extracted from the general 
asymptotic results concerning the oscillatory Riemann-Hilbert problem associated with the second 
Painleve equation (see [IN], [Kapl], [Kit2], [DZ2], and [IK]). 
In the region we define the critical point function as 

*CP {z) = ^ {Rl)-'" V{z) % [N^'Hiz)) , (1.88) 

where the matrix valued function V{z) and the function ^{z) will be defined below in Section 4 [see 
(4.65) and (4.82)]. Both V{z) and C{z) are analytic in 17*^, detV{z) = 1, and C{z) has no critical 
points in Q^. It is worth noticing that the functions V{z), ^{z), and (^{z) depend on the parameter 
y. We will prove that in 

*„(z) = (1 + 0(A^-^))^CP(^), zen"". (1.89) 

We can now formulate the main result concerning the double scaling limit for orthogonal polyno- 
mials. 

Theorem 1.2. There exists 1^2 > such that for all ^2 in the interval < 1^2 < '^'^^ '^^^ 
di satisfying inequalities (1.60), the following holds. Let il", Q} , and J7i he the regions defined 
in (1.61-63) (see Fig. 2). Let Tq > be an arbitrary number and the variable y, defined in (1.43), 
satisfies bound (1.44)- Then the recurrence coefficients Rn obey asymptotic formula (1.45), and for 
the vector valued function ^n{z), asymptotic relations (1.64), (l-^^)> (l-'^S), (1.80), and (1.89) 
hold. 

Remark. It is interesting to compare Theorem 1.2 with the results of Deift-Kriecherbauer- 
McLaughlin-Venakides-Zhou [DKMVZl] and Baik-Deift-Johansson [BDJ], which are both based 
on the Riemann-Hilbert approach as well. The general theorem of [DKMVZl] can be applied to 
our problem at the critical point, y = 0, and for this case it gives asymptotics formulae for the 
function ^n{z) with an error term of the order of A^~^/^. In the work of Baik-Deift-Johansson 
[BDJ], asymptotics similar to the ones of Theorem 1.2 (with the error term of the order of N~'^/^) 
have been obtained for the double scaling limit of a Riemann-Hilbert problem on the circle. We 
provide some comments on the relation of our approach and the one of [BDJ] in Appendix F. 

A different type of the double scaling limit in discrete string equation (1.23), which leads to 
the appearance of the Painleve I equation, was discovered in [BK], [DS], and [GM] in connection 
with the matrix model of 2D quantum gravity. This limit, on the level of solutions of the string 
equation, was analysed in [FIK] in the framework of the Riemann-Hilbert isomonodromy approach. 
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Our method is very close to the scheme of [FIK]. An important difference is, however, that we 
construct an approximation of the orthogonal polynomials of the order of N~'^ and in a finite 
neighborhood of the critical point. The original approach of [FIK] gives only estimates of the order 
of iV-V3 and in a neighborhood of the size of N ' . 

Finally, it is interesting to notice that the Hastings-McLcod solution to Painleve II also appears 
in the Tracy- Widom distribution function at the edge of the spectrum (see [TW2]). 

Double scaling limit for correlation functions. From Theorem 1.2 we will derive the following 
results concerning the double scaling limit for correlation functions of the quartic matrix model. 

Theorem 1.3. Let ^{z;y) = ( ^2^^!^^ 1 the solution for n = to system (1-81) and 

relations (1.85-87). Then the following double scaling limit holds: 

J™ 7 , L m-l ^J^rn i "wo , • • • , ; U + CQi/iV"^/^) 

(ciVVs)™' \cN^/-^ cN'-/-^ ) (190) 

= det (<3c(^ii,^ij;y))ij=i,...,„, 

where c = C'(0) > 0; o,iT'd 

<^^iu;y)<^^{v;y) -^\v;y)<^^iu;y) 

Qc{u, v; y) = r . (1.91) 

it{u — V) 

Furthermore, if z is in the bulk of the spectrum, i.e., < < zq, then 

hm , , , ^ -KATrn\Z-\ , . . . , Z H \tc + CnyJS ' \ 

N^oo {p{z)N)^-^ ^""y p{z)N' ' p{z)N' ' J (1.92) 

where 

^ , , sin7r(n — v) 

Qb{u,v) = ^ 1.93 

Tr{u — v) 

the sine kernel. At the edge of the spectrum, 

1 / 111 Ujji ,^ 2/s\ 

where c = w'{zo) > and 

u — V 

the Airy kernel. 

Remark. The sine-kernel in the bulk of the spectrum is established for a general (fixed) V{M) 
by Pastur and Shcherbina [PS] in a very different approach (see also a nonrigorous derivation of the 
sine kernel in the physical work of Brezin and Zee [BZ]). In the Ricmann-Hilbert approach the sine- 
kernel in the bulk of the spectrum is obtained for a general (fixed) V{M) by Deift, Kriecherbauer, 
McLaughlin, Venakides, and Zhou in [DKMVZl] (see also earlier paper [BI] for the quartic noncrit- 
ical case) . The Airy kernel at the edge is established for the Gaussian model by Bowick and Brezin 
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[BB], Forrester [For], Moore [Moo], and Tracy and Widom [TW2], and for the quartic model by 
Bleher and Its [BI]. 

Plan for the rest of the paper. In Section 2 we will give a formal (perturbative) derivation of 
asymptotics (1.45) for It will be justified later in subsequent sections. In Section 3 we will 
discuss a normal form for the system of differential equations on the psi-function at the critical 
point. In Section 4 we will develop a three step construction of the approximate solution in a 
neighborhood of the critical point. At the critical point four turning points coalesce and this 
makes the asymptotic analysis rather nonstandard. We will be looking for the solution in the form 
V{z)^{N^/^C{z)), where $ (z) is a psi-function for the Hastings-McLeod solution to Painleve II, 
and first we will construct ({z) in the zeroth order approximation, then we will construct V{z) in 
the zeroth order approximation, and after that we will correct ({z) to include terms of the first 
order. The basic role in this construction will be played (cf. [Ble2], [Kap3]) by the equation of the 
equality of periods [see (4.64)]. 

In Section 5 we will construct the WKB approximate solution in Q'^ and we will prove that 
it matches the critical point approximate solution up to terms of the order of N~^. In Section 6 
we will construct a turning points approximate solution in the region $7^, and wc will show that it 
matches the WKB approximate solution. In Section 7 we will construct the WKB solution of cosine 
type in the intermediate region and we will prove that it matches the CP and TP solutions. In 
addition, we will prove that it matches the WKB solution on the top and on the bottom of Hi. 

The analysis of Sections 5-7 provides us with an explicit matrix- valued function ^^i^) which 
solves asymptotically, as iV oo, the basic Ricmann-Hilbcrt problem (1.35-38). In Section 8 we 
will prove that the quotient ^'n(2;)[^° (z)]~^, is equal to /-|-0(A^~^(1 + 12;|)~^) and this will conclude 
the proof of Theorem 1.2. We emphasize that we only use differential equation (1.27) to motivate 
our choice of the function ^^{z). The uniform estimate for the difference ^'„(z)[^'yj(2;)]~^ — / is 
proved by means independent of the WKB theory of differential equations. Indeed, the proof of the 
estimate is based on the analysis of the Riemann-Hilbert problem which is solved by the quotient 
^niz)[^niz)]~^ ■ The ideas and techniques used in Section 8 are close to the Deift-Zhou nonlinear 
steepest descent method [DZl], although there is an essential difference as well. We give more 
details on this matter in Appendix F. Finally, in Section 9 we will prove Theorem 1.3. 

For what follows it will be convenient to make the substitution 

*„(^) = *„(z) ( ^2vr*) ) • ^^-^^^ 

This does not change the Lax pair equations (1.34) and this simplifies equation (1.35) to 

*„+(z) = *„_(z)(j 7). (1.97) 

For the sake of brevity we will redenote '^n{z) back to ^'^(-z) and we will account substitution 
(1.96) at the very end. 
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2. Formal Painleve II Asymptotics Near the Critical Point 

String equation (1.24) is supplemented by the initial conditions 

Ro - 0, i?i - ^-NV(z)dz ^^'^^ 

J — oo 

[see (1.12)]. The critical point Ac = |^ is a bifurcation point for Namely, for A < Ac the 
numbers i?„ are attracted to two different branches for odd and even n, so that 

( R{\) if n = 2A; + l, 
lim i?„ = <^ ; ; 2.2 

n,N^oo: n/N^X \ L{X) if U = 2k, 

where 



H,L= -'^^f^ , A<A„. (2.3) 

sec [BI]. Contrariwise, for A > Ac, the numbers Rn are attracted to one branch, 

lim = i?(A), (2.4) 

n,N—>oo: n/N—>X 

where 



ii= -'^^f^ . A>A. (2.5) 



At the critical point, 

R = L = — — — , A = Af. 
25 




Fig 3: Computer integration of the string equation: i?„ versus n. The solution shown corresponds to 
t = -l,g=l, and N = 400. 
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Fig. 3 presents results of numerical integration of string equation (1.24) with the parameters 
t = —1, g = 1, and N = 400. In this case Ac = 1/4, so that ric = 100. The authors thank 
Bobby Ramsey for his help in carrying out the numerical integration. It is worth noticing that the 
numerical integration has been done by minimizing a global variational functional on trajectories 
{Rn, n = 0,1, . . . , uq} which proved to be an efficient method of integration. 

The problem of the double scaling limit is to find the asymptotic behavior of Rn near Ac- First 
we will calculate the double scaling formally (by perturbation theory) and later we will prove it 
rigorously. Let us assume that 



Rn 



R{y) = -^ + N-l^u{y)+N-'^v{y), if n = 2A: + l, 
L{y) = -^-N-Pu{y) + N-^v{y), if n = 2k, 



(2.6) 



where y is determined by the equation 

Ti 

- = Ae + iV-"y, 

with some exponents a,/3, 7 > to be determined. Our assumption here is that u{y) and v{y) are 
smooth functions of y € M. Now we substitute ansatz (2.6) into (1.24), first for odd n and then for 
even, 

K + N-'^y = R{y) [t + g{2L{y) + R{y)) + g iV^"-^ AL(y)] , 
Ac + N-'^y = L{y) [t + g{2R{y) + L{y)) + g iV^"-^ Ai?(y)] , 
where A is defined as 

/ (y - iV°-i) - 2f{y) + f{y + iV°-i) 



A/(y) = 



For our calculations we can replace A by the operator of the second derivative. Subtracting the 
second equation in (2.7) from the first one, we obtain that 

{R -L){t + gR + gL) + g N^'^'^iRL" - LR") = 0. 

Substituting (2.6) and neglecting higher order terms in we obtain that 

Aguv + Ar2«-2+7 tu" = 0. (2.8) 

To get a nontrivial scaling we set 

2a - 2 + 7 = ; (2.9) 

then (2.8) reduces to the equation 

tu" 

v = -Y-- (2.10) 
4gu 

The first equation in (2.7) gives that 

N-^'y = -N-^2tv - N-'^'^gu'^ (2.12) 
(modulo smaller terms). To get a nontrivial scaling we put 

a = 7 = 2/? 
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Combining these relations with (2.9), we obtain that 

« = 7 = | P=l, (2.13) 

so that (2.6) is written as 

i?n = + i-ir+'N-'/My) + N-'/My), y = n'/' - a^) , (2.14) 

and (2.12) reduces to the equation 

y = -2tv-gu^. (2.15) 
Substituting v from (2.10), we obtain that 

u"{y) 



25 u{y) 

which is the Painleve II equation. Prom (2.10), 



gu\y), (2.16) 



t u"{y) y + gu^jy) , . 

To bring (2.16,7) to a standard form of the Painleve II equation we make a rescaling of u,v and y. 
To that end we rewrite (2.14) as 



Rn = +Ci(-l)"+lAr-V3^(y) +c2Ar-2/3^(y)^ y = ^-lj^2/3 



N~Tg) 



(2.18) 



where 



^2sl/3 /2|t|\^/^ 1/ 1 



CO ={777) ' ci={^) ' C2 = ;t{t^) . (2.19) 



Then (2.16) reduces to 

n"=uy + 2u■^ (2.20) 

which is a standard form of the Painleve II equation (in fact, a particular case of the Painleve II). 
As we have already indicated in introduction, ansatz (2.18) had been suggested in physical papers 
by Douglas, Scibcrg, Shenker [DSS], Crnkovic, Moore [CM], and Periwal, Shevitz [PeS]. 
The equations (2.3) and (2.5) give the boundary conditions 



" (2.21) 
[ ^ 0, if y ^ 00. 

This selects a special solution to the Painleve equation, the Hastings-McLeod solution [HM] (cf. 
Section 1). It has the following asymptotics: 

n(y)= Ai(y)(l + 0(e-(2/3).^'^)), „(y) = y + O (e-(4/3)^^'^) , y^oo, 
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where Ai (y) is the Airy function. Of course, the numbers i?„ in (2.18) do not satisfy the Freud 
equation (1.24) exactly, and we have for them the Preud equation with an error term, 

^ = R^{t + gRn-i + gRn + gRn+i) + O . (2.23) 

More precisely, substitution of (2.18) into (1.24) gives that 



Rn {t + gRn-1 + gRn + gRn+i) " = N'^^'^co [v{y) -y- 2u''{y)] 
+ N-\-ir K'(y) - u{y)v{y)] + O (iV^/^^) , 



(2.24) 



which implies (2.23). 

To get an asymptotic expansion for we use (2.18) with 

oo oo 

u{y) =^N-(^/^^^uj{y), v{y) = ^ iV-(2/3bi;,(y). (2.25) 

i=0 3=0 

Then from the terms of the order of N~'^^^ and we obtain Painleve equations (2.20) on no, v^. 
After that from the terms of the order of N~^^^ and N~^^^ we obtain a linear system on ui, vi; 
then from the terms of the order of iV~^ and N~'^/^ we obtain a linear system on U2, V2, etc. To 
illustrate this process assume that g = 1 and t = —2 (general case can be reduced to this one). 
Then 

Rn (t + gRn-i + gRn + gRn+i) " ^ = A^-^/=^2 ^3 (^;o - 2ul - y) + N-\u'^ - nofo) 



+ iV-4/32V3 (^-Auour +v^ + ^(2< - 16no< + 3^;^)^ (2.26) 

/ 2^/3 \ 

- / ^^^^ ^ ^^^^ ^ __ (-^^^// _ ^^^//^ j + . . . 



Equating coefficients on the right to zero we obtain the equations 

vo = y + 2ul, Uq = uoVq, (2.27) 

and 



2^/3 

2l/3 

VoUl + UqVi = -^{uoVq - VqUq), 

etc., from which we subsequently determine the functions Uj{y) and Vj{y). 



(2.28) 
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3. The ^'-Functions for Painleve II 

Our next step is to derive a model (normal form) equation for the matrix differential equation 

K{z)=NAn{z)^n{z) 

at the critical point z = 0. To that end we substitute (2.18) into the matrix elements of An{z), 
rescale z as 

z = CN-'/h, C=(^^y\ (3.1) 
and keep the leading terms. This gives the equation 

*'(s) = A{s)^{s), (3.2) 

with 

^^^) - y _4s2 + (_i)n2u;(y) - v{y)) -(-l)"4ti(y)s J ' ^"^"^^ 

where 

v{y) = y + 2u^{y), w{y)=u'{y), 

and u{y) is the Hastings-McLeod solution of the Painleve II equation (2.20). For another form 
of the model equation see the physical paper by Akemann, Damgaard, Magnea, and Nishigaki 
[ADMN]. 

Turning Points of the Model Equation. The turning points of (3.2) are the zeros of det A{s) 
on the complex plane. Observe that 

det A{s) = 16s^ + 8ys^ +v^{y) - Wiy), (3.4) 

hence the turning points are solutions of the biquadratic equation 

.^ + p + ^^M^i^ = 0. (3.5) 



Denote for this equation, 



p = p{y)^^, Q = Q{y) = , 

D = D{y) = p\y) - 4g(y) = + 



(3.6) 



Proposition 3.1 (i) The discriminant D{y) is positive for all y G M, so that equation (3.5) 
has 4 roots, ±si^2 = ^si^2{y), such that 

-GO < sl{y) < sl{y) < DO. 
(a) For all y G M, •Si(y) < 0. There exists some yo > such that 

sliy) > 0, y < yo; sl{y) < 0, y> yo- 
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(in) The following asymptotics hold: 



^^'2" 4^4(27r)V2yi/2^ +•••' 



(3.7) 



Proof. We have: 

= - 2u'{y) {yu{y) + 2«=^(y) - u"{y)) = -u\y), 

and D{oo) = 0, hence 



oo 

2/ 



D{y) = / u\x)dx > 0, (3i 
which proves (i). To prove (ii) and (iii) observe that by (3.8), 
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(3.9) 



oo. 



,(0) = -^<0. 

Since q'(— oo) = 0, we obtain that 

I f-y I 
1^y) = lj ^ky)dy=^ + 0{\y\-'^), y- 

This gives the first hne in (3.7). Prom (3.8) and (2.22), 

D{y) = ^e-(^/^)^^^^ (l + O , y ^ oo. 

This gives the second hne in (3.7). Thus, (iii) is proved. 

To prove (ii) we will show that q{y) has a unique zero yo on the real axis. Then (ii) follows. 
Since 

v{y) = y + 2u^{y)>Q, y > 0, 

we obtain from (3.9) that q{y) is increasing for y > and it has one zero yo > on the positive 
half-axis (observe that g(oo) = oo). We claim that g(y) has no zeros on the negative half-axis. 
Indeed, 

v"{y) = 4(u'(y))2 + Au"{y)u{y) = ^{u'{y)f + 4.u'{y)v{y), 

hence if v{y) > then v"{y) > 0. From (2.22) we have that v{y) is negative for negative y 
sufficiently large in absolute value. On the other hand, v{0) = 2u^(0) > 0, hence there is yi < 
such that v{yi) = 0. We claim that v{y) can have only one zero. Assume that there are two. Then 
there exists y2 < such that v{y2) = and v{t) > for y2 — (5 < t < y2, for some S > 0. But then 
'y'(y2) < and v"{t) > iov y2 — 6 < t < y2- Take now any y < y2- By the Taylor expansion, 

1 

v{y) = v{y2) + v'{y2)iy - y2) + ^ J (t - y)v"{t) dt, 
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hence v{y) > as long as v{t) > on [y,y2], which imphes that v{y) > for all y < y2 which is 
not true. The contradiction proves the unicity of the zero of v{y). This in turn proves that q{y) 
has no zero for negative y's. Indeed, q(—oo) = q{yo) = hence if qiys) = for some ys < then 
v{y) = 8q'{y) would have at least two zeros which is not true. This finishes the proof of Proposition 
3.1. 

Standard Form of the Model Equation. The matrix (3.3) differs from the standard Flaschka- 
Newell form of the ^-equation for Painleve II [FN] . To bring it to the standard form we make the 
gauge transformation 

WW = (/*(.), t/=(_'. 7) (J ,_V0 = (-. M)".)- 

Then 

^'{s) = A{s)^{s), (3.11) 

A{s) = U ^A{s)U =[^^^^_ ^ . J , v = y^2u\ (3.12) 

Solutions to equation (3.11) are entire functions of the complex variable s. Their behaivor at 

s = oo is govern by six (the infinity is an irregular singular point of (3.11) of Poincare index 3) 
Stokes matrices. A key fact for our analysis is that the Stokes matrices of (3.11), with u, t;, and 
w determined via the Hastings- McLeod Painleve function, are known. More precisely (see [IN], 
[DZ2]), there exist six Stokes' solutions ^j(s) to equation (3.11) such that 



where 



lim *j(s)e<(^/^)*'+H'^3 ^ J if 



-J 

Isl— >oo 



(j - l)vr 

args — 



<\-e. (3.13) 



3 

The Stokes solutions are uniquely determined by (3.13), and they are related as follows: 

*l(.) = *6(^)(j $2(.) = *l(.)(j J), *3(^)=*2(.), 

$4(^) = $3(^)(_\ J), $5(^) = ^4(^)(J J), MS) = ^^{S). 



(3.14) 



Remark . The existence of the Stokes solutions ^j(s) is a fact of the general theory of systems 
of linear ODEs with rational coefficients (see e.g. [Sib]), and it has nothing to do with the particular 
choice of the parameters u, y, and w in (3.11). For any triple u, y, and w the solutions *j(s) exist, 
they arc holomorphic with respect to u, y, and w, and the asymptotics (3.13) is uniform with 
respect to u, y, and w varying in a compact. Also, in general case, the relations (3.14) are replaced 
by the general equation, 

where the Stokes matrices Sj are some transcendental functions of parameters u, y, and w. A 
remarkable fact ([FN], [JMU]; see also [IN] and indeed classical works of R. Garnier [G]) is that the 
Stokes matrices Sj form a complete set of the first integrals for the Painleve equation (2.20). The 
particular choice of Stokes matrices indicated in (3.14) corresponds to a selection of the Hastings- 
McLeod solution of the Painleve II equation (2.20). 

Equations (3.14) lead in turn to three special vector solutions to (3.2). 
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Proposition 3.2. Assume that A{s) is defined as in (3.3). Then for every real y, there exist 
vector solutions = <i(s;y) and ^j{s) = ^j{s;y), j = 1,2, of the equation $'(s) = A{s)^{s) on 
the complex plane such that: 

(i) $(s) is real, i.e., 

= $(s), (3.15) 

and it satisfies the parity equation, 

= (-l)"a35(s). (3.16) 

As \s\ — > oo, 

(°° -pO \ 

^ ^1 e-'((*/^)*'+^*) , if £ < arg s < TT - £, (3.17) 



(is > 0), with rj] = ( ), and 



r 



oo 



II I arg s\ < — — e or | arg s — 7r| < — — e, 



where o,re some matrix valued coefficients, 
(a) ^i(s) has the asymptotics as \s\ — ^ oo, 



^As)-{f:§]e<^'/'^^''+y^), if -l+e<e.rgs<^-e, (3.18) 

i-ir 



s^ ' 3 ° 3 

Vfe=0 / 



("mj $2(5) = ^i(s) and 



1 -i 

1 y ' 



(3.17') 



^2(-.) = ( ^ ^l^' ) ^i{s). (3.19) 



In addition, 

$i{s)-$2{s) = -i${s). (3.20) 

In what follows we will use the matrix valued functions 

$"(5) = (s), $i(s)) , $^(5) = (s), $2(s)) . (3.21) 
As follows from Proposition 3.2 they satisfy the relations, 

^^{s) = ^'^{s)S, S -- 

$"''«(-s) = (-l)"(73$'''"(s)(T3, l«(i) = ^'^{s). 
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(3.22) 



At infinity for all e > the function $"(s) has the asymptotics 



(3.23) 



as s — ^ oo, e<args<7r — £, 

where 

mr 

7 = , mo = I, 

and the rest of the matrix coefficients mj can be found recursively by substituting series (3.23) into 
equation (3.2). The first two coefficients are given by the following equations: 

iD (-l)"u n^-D^ ^^„w + uD 

mi = --a3-^— a„ m, = I + (-1) a^, ^3^4) 

D = — — yv? 1 
where Uj are the Pauli matrices, 

Proof of Proposition 3.2 is given in Appendix A below. 

4. Semiclassical Approximation Near the Critical Point 

In this section, which plays the central role in the whole paper, we will construct a semiclassical 
approximation to the equation 

^'{z) = NAl{z)^{z\ (4.1) 

in a fixed neighborhood of the critical point z = 0. It is worth to underline from the very beginning 
that we will not prove any results concerning properties and asymptotics of solutions to equation 
(4.1). It will be used only as a tool in the construction of an approximate solution to the RH 
problem. 
Denote 

el^t + gRl + gRl^,. (4.2) 

Prom (1.66) we get that 

el = N-^/^c, + N-'ce + O , C5 = (^1^) ^ [v{y) + 2(-l)"t^(y)] , 

C6=(^) [2{-iru{y)v{y)+v'{y)]; 
e°_, = N-'/'cr + N-'cs + O [n-'/') , cr = [v{y) - 2{-irw{y)] , 

^^=(2^) [2(-ir^(yMy)-^'(2/)]- 
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(4.3) 



This gives the matrix elements of A^{z) as 



«n(^) = - [cig{-ir^'N-'/My) +C29N-^/My) 

aUz) = -{R\ 



z, 



(4.4) 



0x1/2 



where the constants ci, C2 are defined in (1.45). In the semiclassical approximation the function 
detA^(z) will be important. We will use the function 



-,2 ^4 



d{z) ^ -^(.2 - Zi) + - Ae) Z' - {csN-'/' + C4iV-V3) , 

as a suitable approximation to det^°(2;) [cf. (1.58)]. Prom (1.67) we have that 

det A'iiz) = d{z)+a^^z^+f3^„; = 0{N-^/^), /j" = 0(7V-2) 



(4.5) 



(4.5') 



Three Step Critical Point Solution. In a neighborhood of 2; = we are looking for a critical 
point approximate solution to equation (4.1) in the following form: 



*cp (z) = V{z)^ (N'^Hiz)) , 



(4.6) 



where V{z) is a gauge matrix-valued function, ^{z) is a matrix-valued solution to the model equa- 
tion 

$'(z) = A{z)'^>{z), (4.7) 

with A{z) defined as in (3.3), and C(z) is an analytic change of variable with C'(0) 7^ 0. We will 
choose ^{z) differently for Imz > and Imz < to secure the necessary multiplicative jump. Here 
we will carry out a general analysis of (4.1) and we will assume that ^{z) is any solution to (4.7). 
Substituting (4.6) into (4.1) we obtain the equation 



V{z) ^C{z)N-^/^A (Ari/3^(z))] V-\z) = Al{z) - N-^V'{z)V-\z). 



(4.8) 



From this equation we determine itcrativcly ((z) and V{z) in three steps: (1) the zeroth order 
approximation for C(^)) (2) the zeroth order approximation for V{z), and (3) the first order ap- 
proximation for C(^)- 

Zeroth Order Approximation for Ciz)- Taking the determinant of the both sides in (4.8) we 
obtain that 



[C'{z)]''N-^/^detA(^N^/^C{z)) = det [A^z) - N-^V {z)V-\z)] . 



(4.9) 



In the zeroth order approximation we will neglect terms of the order of iV ^ . So we drop the term 
N-^V'{z)V-'^{z) on the right, 

[Ciz)]^N-^/' det A (^N'/'az)) = detAliz), (4.10) 
which is an equation on ({z) alone. By (3.4), 

det A (^N^/^C) = /(C) = 16C^(^) + 8N-^/^yC\z) + iV'^/^ [v^{y) - 4w^{y)] . (4.11) 
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We replace det by d{z) [see (4.5)], reducing (4.10) to the equation 

(Crf{C) = d{z). (4.12) 
Now we drop the 0{N~^) terms in / and d and we introduce the functions 

/°(0 = 16C' + 8iV-^/=^yC'- (4.13) 

and 

d°(.) = ^ (Zo^ -,^)+g(^^- A,) z\ (4.14) 
Equation (4.12) reduces then to 

(C')V°(C) = d°(^), (4.15) 

or 

CVM) = VdPiz). (4.16) 

To fix a branch for the square roots, observe that both f^{z) and d^{z) are positive at z = zo/2 
for large N. We will assume that the both square roots are also positive at zo/2. Equation (4.16) 
is separable and it is easy to solve it. The problem is to find an analytic solution. Consider this 
problem more carefully. Let us make the change of variables 

z = CN-^/^s, C = N-^/^a, (4.17) 

where C is the same as in (3.1). We will call {s,a) local coordinates (at the critical point) and 
{z,() global ones. In the local coordinates equation (4.16) reduces to 



c7'vVV) = V^, (4.18) 

where 



/(a) = A^4/-V"(A^"'/V) = 16^4 + 8ya\ 



where cq is defined in (1.43) and 



(4.19) 



214/3 2/3 

^9 = (4-20) 



Using the definition of y in (1.43) we reduce (4.19) to 



/((t) = 16(7^ + 8?/C7^ 

(5°(s) = 16s^ + 8ys2-iV-2/3c9s6. 



(4.21) 



Observe that ip'^ and S'^ differ only by a term of the order of Therefore, we are looking for 

an analytic solution to (4.18) in the form 

a{s) = s + Ar-2/Vi(s) + Ar-^/V2(s) + . . . (4.22) 

However, there is an obstruction to the analytic solution. Namely, the function [a) has a double 
zero at cr = and two simple zeros at cr = iuo, ctq = (— y/2)^/^. Similarly, the function (5°(s) 
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has a double zero at s = and two simple zeros at some points s = itso that arc close to ^ctq, 
So = o"o + 0{N~^^^). Consider a fixed (i.e., independent of N) closed contour F on the complex 
plane that goes around all these zeros. Then both ^^/ f'^{a) and \/S^{s) are analytic on T, and if 
equation (4.18) holds and a{s) is analytic, then necessarily, 

^ ^/(pO{a)da = j ^/¥{s)ds. (4.23) 

This equation is also sufficient for the existence of an analytic solution cr(s). Namely, integrating 
the both sides of (4.18) we obtain the equation 

/ ^/^(T)dT + Cl= j ^/¥{T)dT + C2. (4.24) 

J CTn J sn 



The integrals on the both sides have the same s^/^-singularities at cr = (Jq and s = sq, respectively, 
and to have an analytic solution a{s) at s = Sq we take Ci = C2. Similar considerations hold at 
s = —Sq and the necessary and sufficient condition to have analyticity of a{s) at both s = sq and 
s = — So is that 

/•—Co r — so 

/ ^/^^dT= / ^/50{^dT, (4.25) 

•J ao -J So 

which is equivalent to (4.23). We will call (4.23) the equation of periods. 

From (4.21) it is not difficult to evaluate the periods by perturbation theory in and 
they are certainly different. Therefore we need a parameter to adjust the periods. To that end we 
change the relation between n and y from (1.43) to 

y = Co'N'/'{^-Xc)+aN-'/^ (4.26) 

where a is a parameter. Then (4.21) changes as follows: 



(4.27) 



and we find the value of the parameter a = a{y) from equation (4.23). Let us analyze this procedure 
in more detail. 

When y = we take a = 0. Then both both ip^ and 6^ have a quadruple zero at the origin 
and equation (4.23) holds. For y ^ 0, let us scale the functions (p^{a) and (5°(s) to the ones 

5{s) ^ S\V^ys) = -s' + (1 - aiV-2/3),2 _ iL_^,6^ ^ ^ 

where y/—2y > for y < and Im \J—1y > for y > 0. Then (4.23) reduces to 

/ VW)dT= I \fi{j)dT. (4.29) 
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where s is the close to 1 zero of d{s), 

S{s) = 0, s = 1 + 0(Ar-2/3) 



(4.30) 



To analyze (4.29) we use the implicit function theorem. To that end, introduce, for small |a;i|, \x2\, 
the function 

I{XI,X2)= / ^-s^ + {l-Xi)s^ -X2SUS, (4.31) 

where s(xi,X2) is the close to 1 zero of the function under the radical. We claim that I(xi,X2) is 
an analytic function at xi = X2 = and 



dI{xi,X2) 



dxi 



7^0, 



dI{xi,X2) 



Xl=X2 = 



dX2 



(4.32) 



Xl=X2=0 



To prove the analyticity apply the Cauchy theorem and rewrite I{xi,X2) as 

I{xi,X2) = ^S ^J-s"^ + (1 — xi)s'^ - X2S^ ds, 



(4.31') 



where Tq is a circle on the complex plane of radius 1 centered at s = 1. In this form one can 
differentiate I(xi,X2) with respect to xi,X2 and prove the analyticity at x-\_ = X2 = 0. To prove 
(4.32) differentiate (4.31'), set xi = X2 = 0, and return back to the real integral: 



dI{xi,X2) 



dxi 



X\=X2=0 



Jo 



2V-s^ + s2 



ds ^ 0, 



and similar for ^^-^ . Equation (4.29) is equivalent to 



dx-2 



/(ari,X2) = /(0,0), xi=aA^-2/3, 



Ar-2/3 



X2 



The implicit function theorem ensures the existence of an analytic solution xi = xi{x2) with 
dxi/dx2{0) 7^ 0. This in turn gives an analytic solution a = a{y) ~ cy^, c 7^ 0, to (4.23). 

Let us summarize our calculations. In the zeroth order approximation the change of variable 
function Co(-z) is determined from the initial value problem 



where 



Co\/7^=v/d^, Co(o) = o, 

d°(.) = ^(^o^-.2)+5(^-A.) .^ 



(4.33) 



(4.34) 



and a is determined from (4.23). We consider the solution to (4.33) in the region 
(-r^i) U U = {z : -zo + di <Rez< zo- di, \Imz\ < ^2}, 
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see Fig.2 above. Wc also notice that, because of the uniqueness of the solution of initial value 
problem (4.33), the following symmetry relation talks place: 

Coi-z) = -Co{z), (4.35) 

When z is separated from 0, i.e. \z\ > coq > where loq < di does not depend on A^, (4.34) 
gives that 

V70(i) = Az^ + Ar-2/3yz + 0{N-^/^), 
,/dP{zj = ^/d(^ + 0{N-^/^), 
[cf. (4.5,9)], hence equation (4.33) implies that 

Kiz)Ciiz)+N-'/^yCo{z) = ^ + 0{N-^/% \z\ > u^. 
In what follows we use the function 

lio{z) = V-diz), (4.36) 
rather than ^/d(z), so we rewrite the last equation as 

4iC^(z)Co'(^) + N-^/HyCi^iz) = i^o{z) + 0{N-^'^), z € ((-^^i) U Q° U n^) \ {\z\ < ojo}- (4.37) 

Observe that the function —d{z) is positive for large iV if 2: > zo + ^^i [cf- (4.5,9)] and we take 
the branch of the square root for /io(-z) which is positive for z > zq + di. In (4.37) we continue 
IJ.o{z) analytically going around zq from above. We note that, for sufficiently large A^, //o(^) is 
holomorphic in the domain 



{|z|>c^o}\((-oo,-z^]U[z^,+oo)), 



where z^ denote the zero of d{z) which approaches zq as N ^ 00. In particular, ij,o{z) is holomor- 
phic in the domain ((— Qi) U U Jli) \ {\z\ < uiq}. 

Equation (4.37) can be used to describe, within an error term of the order of N~'^^^, the 
function Co(^) by an elementary explicit formula. Indeed, as long as z is separated from 0, the 
function ij,o{z) admits the asymptotic representation [cf. (4.5,36)], 



z G A(u;o) = {{-^1) U U Oi) \ {\z\ < loq}. 

where ^J Zq — z^ denote the branch of the root which is analytic in C \ ((—00,2:0] U [2:0,00)) and 
positive for —zq < z < zq. This estimate allows us to rewrite (4.37) as 

v2 



<{zXl{z) + N-^l'y^z) = ^-Jzl -z^ + ^L=iV-2/3 + 0{N-^'% z e A(u;o). (4.39) 

^ * \/Zq - z^ 

By integrating the last equation and taking into account the symmetry z —z we arrive to the 
relation 

^Co'(^) + A^"'/Vo(^) =^oo(^)+iV-2/3yL>i(z) + 0(Ar-4/3), A(a;o), (4.40) 
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where we have introduced the notations 



D^{z) := i^^zl-u^du (4.41) 

and 

D,{z):=co r^i=^- (4.42) 

Notice that both the functions Doo{z) and Di{z) arc analytic and odd in C \ ((— cx), zq] U [zq, oo)). 
(Of course, £)oo(2;) and Di{z) can be expressed in terms of elementary functions, but we will 
not need these expressions; the integral representations are already elementary enough and quite 
convinient for any further analysis.) Equation (4.40), in its turn, implies the estimate 

Co(^) = Coo{z) + iV-2/3yCi(z) + 0(iV-4/3), z G A(cc;o), (4.43) 

where the functions C,oo {z) and (z) are defined by the equations 



Coo(^) 

and 



^^oo(^) 



1/3 

(4.44) 



Di{z) - Coo{z) 

respectively. By a straightforward calculation one can sec that both Coo (2) and C,\{z) are analytic 
at z = 0. In fact, the first terms of the relevant Taylor series are 



Coo(^)=C-iz-^^z3 + ..., Dr{z) = C-^z--^^z' + ..., (4.46) 



and 

(note the absence of the singularity of Ci(^) at 2; = 0). Therefore, the functions Coo (2) and Qi{z) 
are analytic in the full rectangle (— rii) U Jl'' U ili. This, in virtue of maximum principle, allows us 
to extend the asymptotic formulae (4.43) to the full rectangle, 

Co(^) =Coo(^)+iV-2/Vi(^) + 0(Ar-V=^), (-fii)UJ^OuOi. (4.48) 

Initial value problem (4.33) which we have used to define the function Co (-2) was obtained 
by neglecting terms of the order of in basic equation (4.9). Therefore, in view of uniform 
asymptotics (4.46), it makes perfect sense to redefine C,o{z) as 

Co(2) = Coo(^)+iV-2/3y^l(z), y = Co-^iV2/3(^-Ae), (4.49) 

(note that we dropped the parameter a in the definition of y) . With this new definition we preserve 

all the three basic properties of the function Co (-2) which we will use in the next step and later in 
matching the critical point and WKB asymptotics (see section 5 and Appendix C). These properties 
are: 

(i) analyticity and Co(-^) 7^ in the rectangle (— J^i) U 0° U 
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(ii) estimate (4.37), 

(iii) estimate (4.22) for the function a{s) = N^/^Co{CN-'^/^s), 

(iv) symmetry relation (4.35). 

Our next step will be a construction of the gauge matrix V{z) in the zeroth order approximation. 

Gauge Matrix in the Zeroth Order Approximation. In the zeroth order approximation equation 
(4.8) reduces to 

V{z) [c[,{z)N-'/'A [n'/%{z))] V-\z) = A^iiz), (4.50) 
which can be viewed as a generalized eigenvector problem, because by Step 1 the matrices 



CUz)n-'/'a[n'/'Co{z)) 



and A^{z) have the same (up to terms of the order of N~*/^) eigenvalues. We will be looking for 
an approximate solution V{z) to (4.50), with a possible error term in the equation of the order of 
0{N~^), and we will replace the matrix elements of A^{z) by their suitable approximations [cf. 
(4.4)]. We put 

= -aUz) = - \cigi-ir+^N-^/My)+C29N-^/My) 



z . 



aUz) = iRy^'{gz' + N-y'c,) , (4.51) 

Lemma 4.1 Let B = (bij) and D = (dij) be two 2x2 matrices such that 

tr 5 = tr D = 0, det S = det D 

Then the equation VB = DV has the following two explicit solutions: 

Proof of Lemma 4.1 is given in Appendix B below. We apply Lemma 4.1 to solve equation 
(4.50), with 

B = Co{z)N-^/'A [N^/Hoiz]) , D = A^z) 

The problem is that in (4.50) we need an analytic matrix valued function V{z) which is invertible 
in some fixed neighborhood of the origin. Neither Vi nor V2 in (4.52) are invertible. Nevertheless, 
we will find a linear combination of Vi and V2 (plus some negligibly small terms) which is analytic 
and invertible. 

Let us rewrite (4.50) in the local coordinates {s,a) defined in (4.17): 

Vo(s) [a'{s)Aia{s))] V,-\s) = CN^/' Al{CN-"h), 
Vo{s) = l/(CiV-2/3s). 

Define 

B,{s)=a'{s)A{a{s))={^^'^^ ) 

hii = -622 = (-1)"4'U(7'(s)(t(s), (4.54) 
hi2 = (7'(s)(4(t2(s) + (-1)"2u; + v), 
621 = a'{s){-Aa^{s) + (-l)"2u; - v), 
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(4.53) 



and 



dii = -d22 = (-l)"4«s - N-^/^Co\vs + 4s^) 

R'n^9V^\.^2 , . ....... , (4.55) 



To solve (within the error 0{N~'^/^)) the equation 

Vo{s)Bo{s)V„-Hs) = Do{s) (4.56) 

we can use any linear combination of Vi and V2 in (4.52). To determine an appropriate linear 
combination consider the matrix element 611 — dn which appears in Vi and with minus sign in V2- 
By (4.54,55), 

- dii = {-iy'Au(T'{s)ais) - (-l)"4us + N-^/^Co^s{As^ + v). (4.57) 

By (4.22), a{s) = s + N-'^/^ai{s) + ..., hence 

611 - dn = N-^/^Co^s{4s^ + v) + 0{N-^/^). (4.58) 

Thus, the function 6ii(s) — (iii(s) has three zeros: s = and two other zeros which are determined, 
in the zeroth order approximation in N~^/^, by the quadratic equation 

4:s^+v = 0. (4.59) 

(The function 611 (s) — (iii(s) may have more zeros at the distance of the order of N^^^ from the 
origin but we are not interested in those now.) Let us take 

y/detWois) \hiis) - d^^is) 6i2(s)-d2i(s)y ' ^ ^ 

which is a linear combination of Vi and V2 in (4.52). We want Vo(s) to be an analytic invertible 
matrix- valued function. To that end we will slightly correct Wo{s). Consider the matrix elements 
Wij{s) of Wo(s). Prom (4.54,55) and (1.66), 

wii{s) = di2is) - fe2i(s) 

= Ss' + 2v + iV-i/3c-i(-l)"+iu(4s2 + (-l)"2u; + v) + O (^iY-^/^j ; 
w,2{s) = W2i{s) = fen(s) - dii(s) = N-'/^c^'s{4s^ + v) + 0{N-^/'y, (4.61) 

W22{s) = 612(5) - d2l(s) 

= + 2v + Ar-i/3c-^(-l)"+iu(4s2 - (-l)"2w + v) + O (iV-^/^) . 

Observe that in the leading order of approximation in AT"^/^ the zeros of Wij{s) are determined by 
the same quadratic equation (4.59). Let us define now 

~ Wis) - dn{s) - N-y'a2,s b^s) - d2,{s) - Ar-V3«,, ) (4.62) 
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where the numbers an, ai2 = 021, and a-i-i do not depend on s and are chosen in such a way that 
aU the elements in the matrix on the right vanishes at the zeros of 4s^ + v. Then the matrix valued 
function 



8s2 + 2i;' 



is analytic in s and for finite s, 



This implies that 



1 
1 



(2co) 



(2co)-is 
(2co)-^(-ir+^^x 



^ + O (n-'^I^^ . (4.63) 



Fo(s) = 



VdetVFo(s) 
is analytic in s as well and for finite s, 

1 



Wo(s) 



detWo(s) 



:H^o(s) 



K)(s) 



If we go back to the global coordinates z and Co we obtain that 

1 



V{z) = VoiC-^N^/^z) 



^detW{z] 



W{z), W{z) = C-^N-^/^Wo{C-^N^^^z). (4.65) 



To ensure that V{z) is analytic in the region (— f^i) U fi'^ U Oi we will prove the following lemma. 
Lemma 4.2. Define 



W{z) = C-^N-^/^WoiC-^N^/'^z) = 



Wiz) 



Then there exists ci^ > such that if < d2 < then for sufficiently large N for all z G 
(-J7i)Uf^°U0i, deiW{z) 7^0. 

Proof of Lemma 4.2 is given in Appendix E below. Observe that for \z\ < £0 



^^^^ \ aiiiz) - a'ii{z) ai2{z) - a'^iiz) 



(4.66) 



An important feature here is that the correcting terms A'' ^/^an, A 2/'^ai2S, etc., in (4.62) become 
O(A-i) after multiplication by A^-^/s and substitution s = C'^N^/^z. From (3.3), 



Sii(^) = -a22{z) = N-^/\-ir4uCo{z)Co{z), 
ai2iz) = Coiz) [Kliz) + Ar-2/3 ((_i)n2^ + v) 

a2i{z) = Coiz) \-^Coiz) + A^-'/' ((-l)"2w; - v 



(4.67) 



First Order Approximation for Ciz)- We will solve equation (4.9) with an error term of the 
order of A"^. Using notation (4.11) we write (4.9) as 



(C')V(C) = detK(^) - N-^V'{z)V-\z)]. 
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(4.68) 



As a suitable approximation to the determinant on the right we will consider the function 

aiz) = d{z) - iV-i [aUz)q22{z) + 42{z)qii{z) - a\^{z)q2i{z) - a°i(z)(/i2(^)] , (4.69) 
where Qij {z) are the matrix elements of the matrix 

Q{z)^V'{z)V-\z), (4.70) 
where V{z) is defined in (4.65). We replace (4.68) by the equation 

c'/no = xA(i). (4.71) 

In the local coordinates (s,cr) it is written as 

<r'^k{^) = \f^)- (4.72) 

where 

/o((t) = N^'^f {n-''/^u^ = det A{a) = 16(7^ + Sya^ + [v''{y) - 4:w''{y)] (4.73) 

and 

ao{s) = C'^N^/^a (CN-^/^^ . (4.74) 

Prom (4.69), 

ao{s) = do (s) - [dii {s)q^2 («) + d22 (s)g?i (s) - c^i2 (5)^21 («) - c^2i (s)g?2 is)] , (4.75) 

where 



+ [v\y) - Aw\y)\ - N-^'^c^\-ir2w{y) + O [n-^'^^ 



(4.76) 



[use (4.5)], dij{s) arc the matrix elements of the matrix Dq{s) defined in (4.55), and q^j{s) are the 
matrix elements of the matrix 

Qo{s) ^ CN^'^Q (cN-'/h) = V:,{s)V-\s) = 7V-V3 ^^^0 _^ (2^0)"^ ^ + ^ ^^-2/3^ 

[use (4.50)]. Prom (4.75) and (4.55), 

ao{s) = do{s) + iV-i/3(2co)-i [^12(5) + d2i{s)] + 0{N-^/^) 
= do{s) + iV-V3c-i(-l)-2u;(y) + O (iV-^/s) , 

hence by (4.76), 

ao{s) = 16s^+8c^'N^/^ - A,) [v\y) - 4w\y)] +N-^/^rN{s), rjv(s) = 0(1). (4.77) 

Here the notation rN{s) = 0{1) has the following meaning: for any R > there exist A^o = -Ao(-K) 
and C = C{R) such that for all N > Nq, the function rAr(s) is analytic in the disk \s\ < R and 
|''Ar(s)| < C{R) for all s in the disk. 

37 



To get an analytic solution to (4.72) we have to secure the equality of periods of /o(s) and ao{s). 
The function /o(s) has four zeros isj, j = 1,2, (see Proposition 3.1). The function aQ{s) is an 
0(A''~2/^)-perturbation of /o(s) and ao(s) also has four zeros ±tj, j = 1, 2, such that \tj — Sj\ ^ 
as N ^ oo. We have to secure the equality of the periods, 

r ^Us)ds= r ^f^ds, j = l,2. (4.78) 
JO Jo 

To that end we need two parameters. As before we can use the parameter a in (4.26) which slightly 
changes the relation between n and y. Where can we take the second parameter? The idea is to 
change ao(s) by a constant /SN'"^/^ putting 

ao{s) = 16s^ + 8(y- aN-^/^^ + [v'^{y) - Av?{y)\ + N'^/^rNis) + PN'^^^. (4.79) 

Observe that by (4.74), 

therefore, when we will go back to the global coordinates, the function a{z) will change by f3N~'^ 
which is of the order of the error term. We find the values of the parameters a and (3 from condition 
(4.78). Let us analyze this condition more carefully. 

Assume that y yoi so that 7^ 0. Introduce the auxiliary functions 

7j(xi,X2,X3) = / ' Vl6s4 + 8ys2 + \v\y) - 4u;2(y)] + x^s^ + ^2 + x^tn^s) ds, (4.80) 

where tj is the corresponding zero of the function under the radical, j = 1, 2. Then (4.78) reads 

Ij (^-8aiV-2/3,/3Ar-2/3^jV-2/3j =1^.(0,0,0), j = l,2. (4.81) 
Observe that the Jacobian 

J = det 



/ dh 


dh \ 


dxi 


dx2 1 


dh 


dh 1 


\ dxi 


dX2 1 



Xi=X2=X3=0 

is not zero. Otherwise, this would mean that there exists a number c such that the periods of the 
second kind abelian (elliptic in fact) integral, 

I(s) = / — du, 
Jo ^/l6u^ + 8yu^ + [v^y) - 4.w^y)] 

on the genus 1 algebraic curve w = 16s^ + 8ys^ + [v'^{y) — 4w^{y)] are all zero. This in turn would 
mean that the sum, 

makes a meromorphic function on the curve {w, s) having a simple pole at exactly one point, 
namely at one of the two points lying over s = oc. This is imposible since non of these points is 
a Weierstrass point of the elliptic curve (wjS) (see e.g. [FK]). (In fact using Riemann's bilinear 
relations for the periods of abelian integrals, one can show that J = j^.) Therefore, the implicit 
function theorem gives the solvability of equations (4.78) for a and 
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If y = yO) so that S2 = and the curve degenerates, the parameter /? should be taken 

equal to —rf^{0) to ensure that ^2 = as well. The parameter a is determined from equation (4.78), 
J = 1, which is the only period equation left. Its solvability follows from the obvious inequality (we 
recall that yo > 0), 

Jo ^16^^+8^ V 2 

(One can also argue that since the Jacobian J = ^ and hence does not depends on y, the solvability 
for y ^ Uo implies as well the solvability for y = yo-) 

Let us summarize our calculations. The function (^(z) is the unique analytic solution of the 
initial value problem 

C'y/m = V^^ C(0) = 0; 

/(C) = 16C^(z) + 8N-'/^je{z) + [„2(y) _ 4^2(,^)] 

a{z) = d{z) - N-^ [a?i(z)(?22(2) + a'^2{z)qii{z) - al2{z)q2i{z) (4.82) 
-aOi(z)(7i2(z)] -/3iV-^ 
y = Co-^ivV3 (^-A.)+«iV-2/^ 

where d{z) is defined in (4.5), the functions a^^iz) are defined in (4.51), the functions qjk{z) are 
the matrix entries of the matrix Q{z) defined in (4.70) and the constants a and /3 are uniquely 
determined by equation of periods (4.78). The branches of square roots in (4.82) are determined 
for large N by the condition that they are positive on the interval [(2:0/4), (3zo/4)]. The domain 
for ({z) is the rectangle (-O1) U U Oi. 

By exactly the same arguments as in the case of initial problem (4.33) and taking into account 
that the term in the brackets in (4.82) is of the order of iV"^/^ [cf. equation (C.IO) in Appendix 
C], we derive from (4.82) the estimate 

az) = Coo{z) + N-^/^yC^{z) + 0{N-^/^), z e (-J^i) U 0° U Hi, (4.83) 

that is [cf. (4.48)] 

C{z) = Coiz) + O (^N-*/^y z G (-J^i)uQ°UQi. (4.84) 
In addition, the following symmetry relations take place: 

C(-^) = -C(^), (4.85) 

and 

V{-z) = a3V{z)a3. (4.86) 

Equation (4.86) follows from (4.35) and the explicit formulae (4.62), (4.65) defining V{z). The 
uniqueness of a solution to initial value problem (4.82) implies equation (4.85). 

Remark . Unlike the zcroth order change of variable function Co{z), the first-order function 
Q{z) can not be replaced by the first two terms in the right hand side of (4.83). As we will see 
in Section 5 (Proposition 5.2 and Theorem 5.4), in order to match the critical point solution and 
the WKB solution within an 0{N~^) error, the approximation given by (4.83) is not enough. It 
can be, however, used for a more explicit but less accurate approximation to the solution ^n{z) of 
Riemann-Hilbert problem (1.35-38). We discuss this matter in more detail in Appendix F below. 
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5. Matching CP and WKB Solutions 

WKB Solution. The WKB solution is an approximate solution to (4.1) in and it is defined 
as (cf. [BI]) 

^wkb{z) = (5^,r(z)e-['^/>(")''""^+i"~^'"'^^"(")^'(")''"+^^"^] (5.1) 
where Cq ^ 0, Ci are some constants (parameters of the solution), 

^^{z) = [-d{z)]'/\ T{z) = , (5.2) 

and diag A means the diagonal part of the matrix A. Recall that d{z) is a suitable approximation of 
det^°(z) [see (4.5)] and it is a polynomial of the sixth degree in z. Observe that both the function 

1 /2 

lio{z) in (4.34') and ^{z) in (5.2) are defined as [— ii(2)] . The difference is in their domains: for 
11q{z) the domain is ((— Oi) U fi*^ U VLi) \ {\z\ < ujq} while for fi{z) it is Q'^. The function fJ,o{z) is 
an analytic continuation of /u(z) along the contour which goes by the positive half-axis from oo to 
zq + di and then around zq from above. For z lying in [((— Oi) U fi^ U Oi) \ {\z\ < ujq}] H 

f Ho{z), lmz>0, 
I — iJ,o{z), imz < 0. 

The integral /x(u) du in (5.1) diverges at infinity and its regularization is defined as follows. We 
obtain from (4.5) that 

oo 

fl{z) = fXsZ^ + jliZ + /Lt-lZ"^ + Jl{z), Jl{z) = ^ H-3-2jZ~^~'^^ , 

i=o 

q t n 



and we define 



/ u(u) du = —z^ + —z^ + Li-ilnz + / u(u)du. (5.3) 

Because of the logarithmic term, this defines ^("u) du as a multivalued function on the complex 

plane but the function e -'^X^ jg one- valued if n is integer. 

The integral diagr~^('u)r'('u) in (5.1) converges at infinity, because by (4.4,5), as z ^ 

oo, 

[i{z)-a\^{z) z ix{z)-a\^{z) z 

Formula (5.1) can be brought into the following form: 

^'wkb(-^) = Coro(z)e-[^«(^)+^(^)+^^]'^^ , (5.4) 

where Cq = Cq/ \/2 and 

' ai2(«)«2i'(^) - o!i2 {u)a!i^{u) 



C{z) = / fJ'iu) du, t{z) = / 

J oo J oo 



Mu) Hu)-a1^{u)] 



■ du, 
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(5.5) 



The square root branch in (5.5) is determined by the condition that for large z > 0, 

> 0. 



The integral for t{z) in (5.5) converges at infinity, so no regularization is needed. The reduction of 
(5.1) to (5.4) is based on the formula 



diagT-\z)T'{z)=T'iz)as-^ 



In 



/7,(z) - al^jz) 



(5.5') 



which follows from (5.2) by a straightforward calculation (cf. Appendix A in [BI]). The numbers 
Co / 0, Ci in (5.4) are free constants which will be chosen later. 

Proposition 5.1. There exists Nq = No{di,d2) such that for all N > Nq, (5.4,5) define 
^wkb{z) as an analytic function in ff^. As z ^ oo, 

^wkb(^) = V2Co I^I + z-\Rlf/^ (^J J) +0(l^r')) e-[^(^(^)/2)-"•°^+^^l-^ (5.6) 

In addition, the symmetry relation, 

^'wkb(--2) = {-lT(T3^^KB{z)ai, (5.7) 

holds. 

Proof. From (5.3), (5.4) 

*wkb(^) = Coro(z)e-[^(^(^)/^)-"'"^+^^(^)+^(^)+^^]'^% (5.4') 

where ^ 

^(^) = / iW(^) du. 

J oo 

Let us prove that for large N, 

/x(z)-a?i(z) 7^0, zen^. (5.8) 

Indeed, by (4.4,5), 

i.^{z) - {al,{z)f = -d{z) - {aUz)f = ~ - ^ 

+ O (iV-V3(l + 1^1)4) = ^ Q (iV-V3(l + 1^1)4) ^0, Z e 

hence (5.8) holds. In addition, as z — oo, 

\i{z) - a?i(z) 



li{z) 

1/2 



= 2 + 0(|z|-2). (5.9) 



Therefore, the function (^ ^*'^'*^(""^^'* ^ is analytic in il'^. Since fiiz) is analytic in $7"^ and li{z) = 

0{z~^) as z ^ oo it follows that ^(z) is analytic in fi'^. Similar arguments prove the analyticity of 
t{z) and hence the analyticity of *wkb(-2)- Asymptotics (5.6) follows from (5.4' ,9). The equations 

n{-z) = -n{z), a?i(-z) = -a?i, a^i-z) = a^iz), a^i{-z) = a'^i{z), 
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and 

Tl 

imply symmetry (5.7). Proposition 5.1 is proved. 

Critical Point Solution. We define the critical point solution in the region flc = (— r2i)Uf2°Uf2i 

as 

' CV{z)<^'' (n^/^C{z)) , lmz>0, 



. , (5.10) 

CV{z)^'^ In^/^C{z)\ , Im2<0, 



where C is a constant, a parameter of the solution. The functions C(z) and V{z) are defined by 
(4.68) and (4.51), respectively, and the model solutions <I>"'''(z) are defined and described in Section 
3 [see (3.21)]. It is important to notice that equations (4.85), (4.86), and (3.22) yield 

^-cp i-z) = (-l)'^c73^cp {z)a3, (5.11) 

i.e. the same symmetry relation as for the function ^wkb{z). 

Matching CP and WKB Solutions. Let be the horizontal sides of the rectangle Q,c, 

Oiu goal is to show that we can choose constants Co 7^ and Ci in (5.1) and C 7^ in (5.10) such 
that the CP solution ^'cp {z) coincides, up to terms of the order of iV~^, with the WKB solution 
^wkb(-2) on r^. Because of symmetries (5.7,11), it is enough to consider r+. Replacing in (5.10) 
the model function $"'(z) by its asymptotics (3.23), we obtain that 

*Cp(z)=CF(z)(^_^. -^^yo(^)e-<^(4/3X3(.)+;vV3,c(.)-^^)a3^ ^gr+ (5.12) 

where 

+ ^^ + ^^ + 0(iV-^ (5.13) 

To transform ^'cp {z) to the WKB solution we use the following proposition. Denote, as before, by 
Zq the zero of the potential U^{z) (see(1.68)) that approaches zq as ^ 00. 
Proposition 5.2 For z G r+, 



N{A/?,)C,\z) + N^'^'yaz) + N-^/^^C\z)^ =N J\{u) du + t{z) + C° 



+ ^ + iV-2/3A(z)+0(iV-i) 



(5.14) 



where D = D{y) is given in (3.6), 

r°° 1 

C^ = N ^\u)du- -InRl, (5.15) 
/x'=(z) = ^JTP{z) [see (1.68)], and 



^n„.,^/,.^«?l(^) 



^(^) ^ .„0..^ „0.„^ " 2z(-l)"<(z)- 



ia\-^{z) - a%^{z) - in{z) b{z) 
= ' ' i-ir^, b{z) ^ -2pi\z) - ifi{z) {a'i.iz) - a',,{z)) . 
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(5.15') 



In addition, 



where 



V{z) 



1 -i 
-i 1 



To{z)Vo{z) 



(5.16) 



(5.16') 



m = {-ir-a„ n2{z) = -I+ {-iT^'^a^ + ^{z)e{z)a,. 



Proof of Proposition 5.2 is given in Appendix C below. Applying this proposition to (5.12) we 
obtain that 

^cp {z) = c'^o(z)yl(z)e-[^«(^)+"(^)+^°]'^^ 



where 



Yi{z)=Vo{z)Yo{z)e 



Proposition 5.3. Yi{z) = I + 0(iV-^) for z E T+ . 
Proof. From (5. 18, 16', 13) we obtain that 



(5.17) 
(5.18) 



Yiiz) 



I H H 

. az) e{z) 



X 7 + Ar-i/3 



iD 



-as 



_ /v-2/3 

^-l/3ya) ^-2/3y^(2)^ 



az) 

I + A(z)(73 



+ 0(7V-i) 



C(^) 



+ 



+ 0{N-'), 



where 



Y^^^ = m + mi + ^0-3, 



iD 



Y^^\z) = n2{z) + m2iz) + nimi - —I - /\{z)C{z)as + 



(ni + mi)iD 



-(73. 



By formulae (3.24) and (5.16') 

-lYu 



iD 



-cTi, mi 



-era 



(Ti, D = — «^ — yu^. 



n2{z) = -g-^ - (-1)" J + A{z)C\z)as, = / + (-1)" ^ 



nimi = — - 



{-l)"-uD {ni+mi)iD D 



az = —I. 



This gives Y^^"^ = 0, Y^'^\z) = 0. Proposition 5.3 is proved. 

Now we can formulate the main result of this section about the match of {z) with *wkb(-2) 

on r+. 

Theorem 5.4. // we take 



Co = C- Ci = N ^JU^{u)du - - Ini^o, 



(5.19) 
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where the potential function U^{z) is defined in (1.68), then 

*cp {z) = (/ + 0{N-^)) ^-wkbI^) (5.20) 

uniformly with respect to z £ r+. 

Proof. Prom (5.17) and Proposition 5.3, 

*CP {z) = {I + 0{N-^)) C'ro(z)yi(z)e-[^«(^)+^(^)+^°]'^^ (5.21) 

Therefore, if we take the constants Cq and Ci in formula (5.4) as in (5.19), then equation (5.20) 
foUows. Theorem 5.4 is proved. 

Remark. Since the parameters di, d2 of the rectangle CI can vary, estimate (5.20) is valied 
uniformly in an £-neighborhood of P^ . 

6. Matching TP and WKB Solutions 

Turning Point Solution. Let fio = ^^i U (sec Fig.2) and let Cl^''^ be the upper and lower 
halves of Q,o, 

^^''^ = {zeno : ±lmz > 0}. 
The turning point solution *tp (z) is defined in CIq as 



CrW{z)Y^iwiz)), z £ n^, 
CiW{z)Yaiw{z)), z G Q^. 



The constant (7i 7^ is a parameter of the solution. The function of change of variable w{z) and the 
gauge matrix W{z) have been defined by formulae (1.78,79) above. The matrix valued functions 
Y^u,d{z) are model solutions defined in terms of the Airy function, 

Y (.)-(^'^' \(yo{N'/'z) y,AN'/'z)\ 

yuA^)-\^ N-'/^ J \y'o{N^/^z) vl^iN^/^z))' ^^"^^ 

where 

yo{z) = Ai(z), 

yi(z) = e--/6Ai (e-2-/3^), (6.3) 

y2(z) = e'^*/^Ai (e^'^^/^^z) . 

Let us remind that Ai (z) is a solution to the Airy equation y" = zy, which has the following 
asymptotics as 2; — 00: 

Ai (2:) = ^ ^ , ,, exp +0{\z\-^/^)] , -7T + e<avgz<7T-e. (6.4) 



2v^zV4 

The functions yj (z) satisfy the relation 

yi{z) - y2{z) = -iyo{z). (6.5) 

Let us explain (6.1). 
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The turning point solution can be obtained by solving the Schrodinger equation [see (1.50)], 



-ry" + N^U'^ri = 0, 
near the turning point. Namely, we are looking for solutions in the form 



(^N^/Mz)), j = 0,l,2. 



(cf. [Blel]). Equation (6.6) reduces then to the following equation on w{z): 

{w'{z)fw{z) = U'{z) + ^{w,z} 

where {w, z} is the Schwarzian, 

{w,z} 

Dropping the Schwarzian term we get 



w"'{z) 3 fw"{z) 



or taking the square root, 



w'{z) 2 \ w'{z) 
{w'{z)fw{z) = U'^{z), 

l^lw'/\z))' = Vim 



(6.6) 
(6.7) 

(6.8) 

(6.9) 

(6.10) 
(6.11) 



To secure the analyticity of w{z) &t z = Zq we take a solution as 



w{z) 



2/3 



(6.12) 



The function U^(z) is positive for z > Zq and we take the positive branch for ^^/U^{u) for u > Zq 
and also the positive branch for the power 2/3. Formula (6.12) defines w{z) as an analytic function 
in ilo for large N. Observe that because the dropped Schwarzian term in (6.8) is of the order of 
(6.7,12) solves (6.6) with an error 0{N-^). 
To obtain the gauge matrix W{z), recall that equation (6.6) was derived from the system 



^1 

•02 



by solving tp2 in terms of ipi , 



i^2 = — {N Vi - oiiV'i) 



'■12 



and by substitution ipi = \/'A.2 V [cf- (1-49)]. For ipi and ^2 we obtain from (6.7) the following 
formulae: 



V'i(z) =CiN^/^ 



w'{z) 



w'{z) 

1/2 



1/2 



(iV2/3u;(z)) , j = 0,l,2; 



(6.13) 
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(we omit the term of the order of N in -02), or in the vector form, 
where 

W{z)=(^Y'( aii.) (6.15) 

\w'iz)j y-^o^) 

This gives the gauge matrix W{z). 

It is important to notice that the turning point solution *tp {z) has the "right" jump on the 
real axis: 

*+p(z) = *Tp('2)(j 7)' ^eOon{Imz = 0}, (6.16) 

which follows from equation (6.5). We want to check that if we choose appropriately the constants 
Co and Ci in (5.1), then the turning point solution *tp {z) matches the WKB-solution ^'wkb('Z) 
on the boundary of $7o; excluding a neighborhood of the interval [0, zq]. Let Tq be the union of the 
three sides of the square excluding the one crossing [0, zq], 

Fq = {z G r^o '■ llm-^ol = <^2 or Rez = + <^i}- (6-17) 
Lemma 6.1. If we take the WKB solution with the constants 



Co 



Ci 



27ri/2 ' 

then 



r°° , 1 

Ci = N ^U^{u)du- -luRl, (6.18) 



*TP {z) = {i + 0{N-^)) *wkb(^), z e rs, (6.19) 

uniformly with respect to z eTq. 

A proof of Lemma 6.1 will be given below in Appendix D. It will be based on an alternative 
form of the WKB solution in fi'^. 

Alternative Form of WKB solution. The standard WKB form for a solution of the Schrodinger 
equation (6.6) is 

{x'{z)) 

Equation (6.6) reduces then to the following equation on x{z): 

(x'(z))' = [/°(z) + ^{x,4, (6.21) 

where {x, z} is the Schwarzian derivative. Dropping the term with the Schwarzian derivative we 
obtain 

>,'{z) = {U"{z)Y^\ (6.22) 

which gives 

,(.)= Co <A^i:>(."(.))--wc..) 
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and the WKB solution in the form 

^'wkb(^) = CoT-{z)E{Ne{z)), (6.24) 
where the gauge matrix T'^{z) is defined in (1.70), 

e{z) = r fi'iu) du + C^, ^i'iz) = (C/°(z)) , (6.25) 
and the model solution E{z) is 

^(^) = (_'e- e^)- (6.26) 

Lemma 6.1 is an obvious corollary of the following result which will be proved in Appendix D below. 
Lemma 6.2. // we take Co, Ci as in Lemma 6.1 and Ci = then for z G O'^, 

^'wkb(^) = (1 + 0{N-^\z\-^)) ^'wkb(^), z e n^, (6.27) 

and for z eTq, 

*TP {z) = {! + O(Ar-i)) z e r^. (6.28) 



It is worth noticing that the derivation of (6.28) will be based on WKB type asymptotics for 
the Airy function. In the vector form the latter is formulated as follows. Let 

(-) = ( ti^rl V m=( C ) > Ciz) = ( ^ ) . (6.29) 



Ai'(z)y ' \-e~^ J ' ^' \-smz 

Proposition 6.3. For any e > 0, as z ^ oo, 

Ai{z) = [l + 0{\z\-^/^))^B{z)E(-^\, -;r + £<argz<7r-£, (6.30) 

where 

B{z)=(^'~'J' (6.31) 

and 



Ai i-z) = (l + 0{\z\-'/')) ^Biz)C 



- / 2^3/2 27r 27r 

— + £<argz< — -e. (6.32) 



3 4/' 3 -°-3 



The notation 0(|2;|~^/^) in (6.31,32) means a 2 x 2 matrix valued function r{z) such that 
\f{z)\ < Co\z\~^/^ for \z\ > Ci, where Co,Ci > are some constants. The branches for fractional 
powers are fixed by the condition that they are positive on the positive half-axis. We will not prove 
Proposition 6.3 because it is just a reformulation of well-known asymptotics for the Airy function. 
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7. Matching CP and TP Solutions 

In this section we will define a WKB solution ^'wkb(-2) in the region f^i and we will show that 
it matches both (z) and (z). We begin with introduction of an auxiliary CP solution. Let 
Q,^ be the part of flc = OP (J J7i between the diagonals of the first and fourth quadrants, 

= {-2^ € f^c '■ \iToaz\ < Rez}. 
We define the auxiliary CP solution in JIJJ as 

*CP (^) = CViz) [$2{N'/\{z)), 4^{N'/\{z))) (7.1) 

where ^i{z) and ^2{z) are the Painleve II psi-functions defined in Proposition 3.2. Using (5.10) 
and (3.21) we can relate '^^p (z) to *cp {z) as 

^cp (z) = ^hp i^)Si,2, ±Imz > 0, (7.2) 

where 

Observe that ^'^p {z) is defined in terms of ^i{z) and ^2iz) for which we know from (3.18) their 
asymptotics in the sector 

-(7r/3) + £ < argz < (tt/S) - e. 

This allows us to find asymptotics of ^Qp{z) in Qi, assuming that d2 is small enough. Namely, 
C'(0) > £ > and if ^2 is small enough then z e 0,i implies that 

-(7r/3) + £ < argC(2) < (tt/S) - e. 

Using the asymptotics of <&i,2(-2) we obtain the following result. 

Proposition 7.1. If we take Co and Ci as in (5.19) then uniformly with respect to z E fli, 

*Sp(z) = (7 + 0(iV-i)XKBW5o, ^o=(o °), (7.4) 

where the function ^wkbI-^) '^'^ analytic continuation of the function "^wKsiz) in (5.4) from the 
upper half-plane to Oi. 

Proof of Proposition 7.1 repeats the one of Theorem 5.4 and we omit it (take into account 
that formula (5.14) is extended to Qi, see Appendix C). Formula (7.4) is easy to check without 
any calculations if Imz > e > 0, z G fli. In this case $i{N^^'^(^(z)) is exponentially small, while 
$2{N^^^C{z)) is exponentially big (cf. Proposition 3.2), hence in (7.1) the function 

$2{N^^\iz)) = $i{N^/^C{z))+i^{N^^\iz)) 

can be replaced by i^{N^^^({z)) with an exponentially small error. This gives \I'cp (^;)<S'o which 
can be further replaced by '^wkb{z)So, with an 0{N~^) error, due to Lemma 6.1. Thus, we obtain 
(7.4). 

We will call the function on the right in (7.4) the auxiliary WKB solution, 

*wkb(^) = *wkb(^)5o, (7.5) 
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so that 

*£p(z) = (/ + 0(iV-i))^'WB(^), 
With the help of ^wkbC-^) introduce the WKB solution in fii by pattern (7.2), 

*wkbW = *wkb('2)'9i,2, ±lmz>0. (7.6) 

It shares the following nice properties: 

Proposition 7.2. Uniformly with respect to z e fli, 

*cp {z) = (1 + 0(Ar-i))*wKB(^). (7.7) 

On the real axis, 

*wkb(^) = *wkb(^)(J di<z<zo-di, (7.8) 

Finally, on the horizontal sides o/Oi, 

*wkb(^) = (1 + o(iv-^))^wKB(^)> ^ e r+ u rr, 

w/iere 

= G i7i : lmz = ±^2}- (7.9) 

Proof. Equation (7.7) follows from Proposition 7.1; (7.8) from (7.6); and (7.9) from (7.7) and 
Theorem 5.4. 

We introduce next the auxiliary TP solution by the same pattern as the CP one, 

*TP (Z) = CiW{z)Ya{w{z)), Z e nl 

where 

/ArV6 \ /2/2(iV'/'^) yiiN^/^'z) 



The function ^'^p (z) is defined in the domain 

^0 = {-^ = ^0 '■ |Iniz| < ^0 — Rez}, 

between two 45°-lines through zq to the left. Using (6.1,2,5) we can relate 5'|^p (z) to \1'tp (z) as 

*TP {z) = *TP {z)Si,2, ±Imz > 0. (7.11) 

By (6.3,4) we know the asymptotics of yi,2{z) in the sector | arg^; — 7r| < 2tt/3 — e and from this 
asymptotics we obtain, like in Lemma 6.1, that ^'^p {z) matches ^wkb('^) ^i' 

M/^P (z) = (/ + 0(7V-1))^'Wb(^), z e 

This, together with (7.11), proves the following addition to Proposition 7.2: 
Addition to Proposition 7.2. Uniformly with respect to z e Oi, 

^-TP (z) = (1 + 0{N-^))^^kb{z). (7.12) 
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8. Proof of the Main Theorem 

Define the matrix- valued function ^'^(z) on the complex plane by the formulae: 



(-1)"(73*wkb(-2;)(T3, z G (-J^i); 

^Tp(-z), zG^i; 
(-l)-a3^'Tp(-^)a3, zg(-J^i); 



Let 



Prom (1.37) and (5.6) we obtain that as z — oo, X^iz) admits the asymptotic expansion 



j=0 



zi 



with 



In particular, 



lim Xn{z) = Go. 

Prom (8.1) we obtain that *° (^;) satisfies the following parity equation: 

<(-^) = (-l)V3<(z)a3, 
the same as '^n{z) [cf. (1.42)]. Hence 

Xn{-z) = a3Xn{z)a3. 



(8.1) 



(8.2) 



(8.3) 



(8.4) 



(8.5) 



(8.6) 



(8.7) 



The function Xn{z) has multiplicative jumps on a number of contours, because of the piecewise 
definition (8.1) and because of the jump of ^n{z) on the real axis. Let us discuss this question 
more carefully. Let 70 be the union of boundaries of the regions il, ili, and —ili, 

Then by Proposition 7.2, 

[X„_(^)]-iX„+(^) = = I + 0{N-'), z e 70. 

The function Xn{z) has no jump on the segment [—zq — di,zo + di] because here the jumps of 
^ni^) "^niz) cancel out each other. Consider also the set 



7i = {z < -zo - di} U{z> zo + di}. 
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Fig 4: The contour 7 = 7o U 71 . 

and let 7 = 70 U 71 (see Fig.4) . 
On 71, 

Observe that ii z > zq + di then 

where the matrix T{z) is uniformly bounded, hence 



[X^.{z)]-'X^+{z) = [f{z)]- 







1 



with some a > 0. Thus, we obtain the following proposition. 

Proposition 8.1. The function Xn{z) is an analytic matrix-valued function on the complex 
plane with multiplicative jumps on the contours 70 and 71 such that 



I + 0{N-^), ze7o, 



(8.8) 



At infinity Xn{z) admits asymptotic expansion (8.3). 

Proposition 8.1 implies that the fucntion Xn{z) solves the Rimann-Hilbert problem on the 
contour 7 = 70 U 71 , 

Xn{oo)= lim X„(z) = eo (8.9) 



Xn+{z) = Xn-{z)G{z), ZG7, 
with the jump matrix given by the equations, 

G(z) = [<_(z)]-i^o + (2;), ZG70, 



(8.10) 



^.111 



and satisfying the estimates. 



\\I-G{z)\\L,^^)nLUi) = 0{N-'). 
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(8.12) 



The Riemann-Hilbert problem shares a remarkable property of wellposedncss (sec e.g. [BDT], 
[CG], [LiS], and [Zh]). Namely, estimate (8.12) imply the following estimate of Xn{z) on the full 
complex plane (see Appendix D in [BI]). 

Proposition 8.2. For all z e C, 



Comparing (8.13) with (8.3) we obtain that 

e^'e, = o{N-^), 

that is, according to (8.4), 

g-(Ci-A„)^3r-irie(^i-^'')'^3 - (i?°)^/Vi = O(iV-i). (8.14) 

This is a 2 X 2 matrix equation. If we write it for the matrix elements we obtain from (1.38) that 
two equations are trivial, = 0, and the other two are 

(8 15) 



Observe that by (4.2), 



so that > is uniformly bounded and separated from as AT — > oo, hence from the first equation 
in (8.15) we obtain that 

gA„-Ci ^ (i20)i/4 ^ o(iV-i). (8.16) 
Substituting of this into the second equation in (8.15) gives that 

Rn = Rl + (iV-1) . (8.17) 

Prom (8.16) we also obtain that 

e^" =e^HO^/^ (l + 0(iV-i)) . 
Since hn = e^'*'" (see Proposition 1.1), we obtain that 

/t„ = e2^H^n)'^'(l + 0(iV-i)). (8.18) 

The constant Ci is given in (5.19). Substituting its value into (8.18) we obtain that 

K = e •''0 (l + 0(iV-i)). (8.19) 

This proves the part of Theorem 1.2 concerning the asymptotics of i?„ and hn- 
Let us take ^ 

Co = 21/2(^0)1/4- (8-20) 

Then from (8.16) and (8.4) we obtain that 

eo = / + o(Ar-i), 

hence 

Due to (8.1) and (1.96) this proves asymptotic relations (1.64), (1.76), (1.80), and (1.89). Theorem 
1.2 is proved. 
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9. Universality 

Critical Point. The double scaling limit at the critical point z = is determined by the kernel 

1 / 11 V \ 

Q^{u, V) = ^lirn^ ^ j , (9.1) 

where 

^ , ^ „l/2i^N{z)i^N-l{w) - 'lljN-l{z)'ipN{w) 

z — w 

and c > is a normalizing constant. To evaluate (9.1) we apply CP solution (1.88) which gives 
that modulo 0(iV-V3) 

terms, 



V'^(^) = ^ «)'^' \^,,{z)^\N''\{z)) + V,2{z)^\N''\{z)) 
^M-i{z) = ^ «)-^/^ [v2iiz)<P\N'/'az)) + V22iz)<^\N'/'C{z)) 
To evaluate (9.1) we can replace Vij{z) by Vij = Vij{0) and ({z) by C'(0)'2^- We take 

c = C'(0). 

Then, modulo 0(iV-V3) 

terms, 
1 / « 



(9.2) 



VciVi/3' cAfi/3, 

= C(U - V)-^ [{Vu^\u) + yi2$'(u)) {V21<^>\v) + V22$'(l')) 
- (yil$'(f) + Vi2<^^{v)) {V2l'^\u) + V22^'(^i))] 



(9.3) 



$(z) = i7oi?o(^), C/o={ S-i / 'Ii 



u — 

where C = 1/tt and 

C' = C(VnF22 - 142't^22) = Cdety = C. 

Observe that as they are defined in Section 3, the functions ^^{z) depend on n mod 4. Neverthe- 
less, the combination 

$^(u)$2(t,) _ <^^{v)^^{u) (9.4) 
does not depend on n. Indeed, by formula (A. 5) in Appendix A below, 

i-i) 

where .60(2;) = ( ^2^^^ ) does not depend on n. A direct computation gives that 

^\u)^\v) - ^\v)^\u) = (det C/o) [B',{u)Bl{v) - B',{v)BI{u)] . 
Since det Uq = 2i, this proves that combination (9.4) does not depend on n. Thus, 

fe(.,,)=-^'W^^W--^>'^^'">. (9.5) 

Tr{u — V) 
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Edge of the spectrum. The double scaling limit at the edge is determined by the kernel 



1 / 1i V \ 



(9.6) 



where Qn{z,w) is the same as in (9.2) and c > is a normahzing constant. To evaluate (9.6) we 
use the TP solution 



i^n{z) 
i^n-liz) 



OU/4 



27ri/2 
27rV2 



(9.7) 



l^ii(z)iVi/6 Ai {N^/^C{z)) + Wi2{z)N-^/^ Ai'(iVi/3^(z)) , 
VF2i(^)iV'/^ Ai (N^^Hiz)) + W22{z)N-^/'' Ai'(Ari/3^(2))" . 

If we take c = C'i^o) repeat the above calculation at the critical point, then we will obtain that 

Ai(n) Ai'(v) - Ai(t>)Ai'(u) 



Qe{u,v) 



u — V 



(9.8) 



the Airy kernel. 

Bulk of the spectrum. The double scaling limit in the bulk is determined by the kernel 



1 / U V \ 

Q{u, V) = Jirn^ —Qj, + ^ + ^ J , 



(9.9) 



where c = p{z), the value of the density function at z. To evaluate (9.9) we use the WKB solution 
and in the same way as above we obtain that 



Qbiu,v) 



sin7r(w — v) 
7r{u — v) ' 



the sine kernel. 
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Appendix A. Proof of Proposition 3.2 

Proof. Let ^j{z), j = 1, 6 be the Stokes' solutions to equation (3.11) introduced in Section 
3. We recall that the solutions ^^jiz) are uniquely determined by the asymptotic condition, 



lim <ifj{z)e<^^^^^''+y^>'=I if 



are z — 



U - l)vr 



TT 



and that they are related as follows: 
^i{z) = ^g{z) ' ^ ~^ 

^4{z) = ^3{z) 



1 
1 

-1 1 y ' 



^5{z) = ^i{z) 



1 1 

1 



*3(2) = ^'2(2), 



(see (3.13) and (3.14)). Let us rewrite these equations in terms of vector solutions. Let 
Then by (A.2), 



(Al) 



iA.2) 



*l(^) 


= ^i{z) 


^2iz) 


= ^3iz) 


*5(2) 


= ^6{Z) 



(^Boiz),B,{z)) 



B2iz) = Boiz) - Bi{z), 



{A.3) 



and 



Define 



Then 



lim ^o(z)e<("/3)^'+J'^ 



lim Bo{z)e-<^^/^>'+^') -- , ^ 

|z|— >oo V 1 



lim B,(z)e-<^^/^>"+y') -- , , 



lim 52(z)e<(^/3)z^+yz) 

|zHoo 



${z) =i''UBo{z) = 
$i{z) = i''-^UBi{z) 

$2{Z) = i"+^C/^2(2) 



if £ < arg z < n — e, 

if TT + £ < arg z < 2Tr — s, 



.„ TT 47r 

if -- + £<argz< y-e, 

if Y + ^ < a^g^ < Y - ^- 



n-l 



,n+l 



i?l(^), 



(-i) 

^-K:l ]B2iz), 



{-if 



$'Jz) = Aiz)$Jz) 



{AA) 



{A.5) 



iA.6) 
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and 



^ lim %{z)e<^^/^^''+y') = ^ , if £ < arg z < TT - £, 

lim $(z)e-^((^/3)^'+^^) = ( , ^T']!i y if TT + £ < arg z < 27r - £, 

Vi-^J / ^^^^ 

^jim^#,(.)e-<(V3).3.,.)^^H^^^^ .f _|+,<arg.<f -£, 

lim $2(^)e^((^/=^)^'+^^) = f y if ^ + £ < arg 2 < ^ - £. 



Observe that the function $ (z) satisfies the equation (A. 6), because the matrix A(z) is real, and, 
in addition, $ (^) has the same asymptotics (A.7) as $(z). Hence (3.15) holds. Similarly, $i {z) 
satisfies the equation (A. 6), and it has the same asymptotics as ^li.^). Hence (iii) holds. 

/(_!)" \ 

The function z), where ^ ~ ( q ^ ) ^^^^^^^^^ (-^-6) and it has the same 

asymptotics as ^{z). Hence (3.16) holds. Similar argument proves (3.19). Finally, (3.20) follows 
from (A. 3). Proposition 3.2 is proved. 

Appendix B. Proof of Lemma 4.1 

Direct computation: 



and 



Since 



di2 \ / 6ii bi2 \ _ ( hidi2 bi2di2 

^11 - i'lic^ii + bi2b2i -bi2dii 





612 ^ 




U21 


-611 J 





dii di2 \ ( di2 \ _ / biidi2 bi2di2 

2 
11 



£^21 -dii ) \bxi- c/11 bi2 ) \d^i- budxx + £^12(^21 -&i2(^ii ) ' 



b\i + 612621 = - det i? = - det I? = d\i + ^12^21 , 
we obtain that ViB = DVi which was stated. Similarly, V2B = DV2. Lemma 4.1 is proved. 
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Appendix C. Proof of Proposition 5.2 

Let z G r+. By (4.37) and (5.2'), 

p.{z) = 4iC^(z)Co'(z) + N-^/HyCo{z) + O , (C.l) 

hence from (4.66,67) we obtain the following asymptotics of the matrix W{z) as N ^ oo: 

+ N-'-"' [-(-i)"2t«c;(2)<T3 + (t. - !/)c;(2)'i + o (jv-i) . 

Substituting this formula into (4.65), we obtain the equation, 
where 

b{z) = -2^?{z) HAz)[aUz) - al^{z)\, {C.2') 
and the matrix- valued function Uo{z) admits the asymptotic representation. 



Uo{z)=I + N-'l\-ir——a, 

Ko{z) 



+ N- 



-2/3 



+ O (AT-i) . 



(C.3) 



When deriving (C.3), we use (C.l), the formula 

b{z) = -2i^{z)[a%{z) - i^i{z)] + O (iV-2/3) , 
which follows from (C.2'), and the formula 

which follows from (4.4). Taking into account that by (4.5'), 

^i\z) = {al,f{z)+aUz)al,{z) + 0{N-^/^), 
a simple algebra shows that 

^)[-aUz)+iaUz)-t^{z) za?,(.)-aO,(.)-zM.) +^^^ ^^"^ ^ 

where 

_ ( ^i2iz) - ¥{z)+iaUz) Y' ^ I «?2(^)+«^i(^) 



m?i (2;) - (z) - i)u(z) y V (^) - «2i (^) - ^K^) 
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and To{z) is the WKB gauge matrix [see (5.5)]. Prom (C.4), 



S{z) = l + N-^/^ — 



+ o(iv-i) . 



- 0^1(2:) - ifi{z) 

Using this equation, equations (C.3) and (C.4'), and also equation (4.84), we can rewrite (C.2) as 

1 - 



Viz) 



-i 1 



To{z)Vo{z), 



(C.5) 



where 



u 



and 



Voiz) = 1 + Ar-V3(_i)"__^, + iV-2/3 
+ 0{N-'), 

ia?! (z) - (z) - ^//(z) 



J + (-l) 



n+l 



4C2(^ 



-(72 + A(2;)(J3 



- 2z(-l)"<(z 



biz) ■ 



iC.6) 



{C.&) 



This proves equation (5.16). Let us prove (5.14). 

By (4.68,69), 

a{z) = det [Al{z) - N-^V {z)V-^ [z)] + O [N-^) 
Using (C.5) and the diagonahzing property of the matrix Tq{z)^ 



we obtain that 



det - N'WV-^] = det - N'^TI^Tq^ - N'^ToV^V^-^T^-^ 

= det [T^'A'iTo - N-X-X " A^"'V^o 
= det [-/x(73 - N-^T^X - ^"'^0 V'] 

= -/x^ det / + -N-^a^ToX + -iV" VaFo'^o"^ 

= -I? - hN-^Tt [asT^X] - I^N-'Tt [^3^0' V'] + O (iV'^) 

Prom (C.6) we conlude at once that 

Tr [aaV^Vo'^] = 2A'{z)N-^/^ + O {N'^) . 

At the same time, equation (5.5') and the identity (following from the equations detTo 

(ro)ii = (10)22), 

diagTo-^T^ = Utv [asT^X])<^s, 



iC.7) 



{C.8) 
2 and 



imply 



Tr [asT^X] = 2r'(2:). 
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iC.9) 



From (C.7-9), and taking into account our convention about the branches of the square roots in 
(4.82), we derive the asymptotic formula, 



[det{A'i{z)-N-'V'{z)V-\z))+0{N-')f' 
= -ifi{z) - iN-^T'{z) - iN-^/^A'{z) + O , 



(C.IO) 



Prom (4.11) it follows that when z G r+. 



1/2 



V7(0 = {l6C" + SN-'/'yC + N-^/' [v\y) - 4w\y)] } 

= 4C^ {i + iV-/3^ + ^-4/3 ^^(^)-4gy)-y- ^ ^ 

D 



(C.ll) 



2C 



where [cf. (3.6)] 



D = D{y) 



y"^ -v'^{y) + 4w^{y) 



Replacing the both sides of (4.82) by their asymptotics (C.10,11) and integrating the resulting 
equation, we arrive to the formula. 



l^le + N-'/'yC + N-'/'^C') = l\{u)du + N-'r{z) 

+ N-^/^A{z) + c + 0{N-^) . 



{C.12) 



Let us find the constant c. We will use the fact that the function (^(z) is odd in z. 

Consider the contour 7 which goes along the imaginary axis from zoo to id2, then around zero 
on the right to —id2 and then to —zoo. Denote by 71 the part of 7 from ioo to id2, by 72 the part 
of 7 from id2 to —id2, and by 73, from —id2 to —00 (see Fig.5). 































^3 





Fig 5: The contour 7. 



Let iJ,o{z) = ^J—d{z) with two cuts and {z^ ,00). We will consider iio{z) outside 

of a small disk D^^^^ = {\z\ < ujq] around the origin. For z G $7^^, 

JM^), lmz>0, 

/^o(^) = { / X r / n V<^-^^) 
I — /x(z), Imz < 0. 

59 



In contrast to fi^z) which is an odd function, the function /Uo(-2) is even. The function t{z) is defined 
in (5.5). Let tq{z) be the function on C \ D^^^ defined by the same formula (5.5) with two cuts 
{—oo,—z^) and {z^,oo). Replacing in all formulae (C.1-11) the functions //(z) and t{z) by the 
functions fio (z) and tq {z) respectively, we can extend the validity of these formulae to the points z 
of the annulus Clc \ -D^o- Hence equation (C.12) can be extended to the equation 



+ Ar-5/3Ao(z)+c+0(Ar-2), 



{C.12') 



where 

Ao{z) = A(z)|^(^)^^„(^), 

and it is an analytic continuation of A{z) to C \ D^^g with two cuts (—00, —z^) and {z^ ,00). 
We will split the function to{z) as 

To{z) = fo{z) + Toiz) + CO + O(iV-i), (C.14) 

where 



fo{z) = ^In 



aUz) 



°° ^ (C7.15) 



U\z) = a',,' (z) - al.iz)^ 



[we will justify (C.14) below, in the end of this appendix]. The branch of the logarithm in fo{z) is 
determined by the condition that for sufficiently large N, 

0<Imfo(z) < |, zeC\[D^,[J{-oo,-z^]U[z^,oo)]. 



(Note that ^"^^jrg;^"^ = (Piz) + 0{N-^/^), where 4>{z) = and : C \ [(-00, -zo] U 

[zo,oo)] HH- {Re(/) > 0}.) Using (C.14) we rewrite (C.12') as 

,(ie.N-v^yc +^-fr') ^£ (..w.iv-.|M),„,^-.,„(,) ^^^^^ 

+ iV-^/=^Ao(z) + ci + O , ci = c + N-^co. 

where z ^Q,c\ D^^^ . The integral l^oiu) du diverges at infinity so a regularization is needed. Let 
us describe the regularization we use. 

Let 7o be any contour which goes from 00 to z in the region 

So = C\ {D^„ U (-00, -z^) U (z^, 00)) , 

and which starts at 00 in the upper half-plane. For our purposes it will be useful to introduce a 
slightly more general regularization than (5.3), which in fact will be equivalent to (5.3). Take any 
point a on 70 and define the regularized integral as 

/ IjLo{u) du = ^a'^ + + fi-ilna + / 'jlo{u)du+ / iio{u)du, (C'-17) 
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where 

Jio{z) = no{z) - fisz^ - l^iz - n-iz~^; = 2' '"1 = 2' '^~^^~N' 

The branch of In a is taken as follows. By our assumption 70 starts at 00 in the upper half-plane. 

If the whole piece of 70 from 00 to a lies in the upper half-plane, then take In a on the main branch 
of logarithm, with a cut at (— oo,0) and Inl = 0. For other a's extend In a continuously along 70. 
Definition (C.17) shares the following important properties: 

1. The right-hand side of (C.17) does not depend on a. 

2. The contour of integration 70 can be deformed by the Cauchy theorem. 

3. The integral l^o{u) du does not depend on the contour of integration in Sq. 

Observe that the property 3 is not automatic because Sq is not simply connected. It follows 
obviously from the fact that ^o(^) is an even function. 

The function has the following asymptotics as z — 00, Imz > 0: 

^ = -i.-'+rW, r(.) = 0(.-»), (C.18) 

[see (1.71)], and we define the regularized integral of this function as 

/ \-^du = h / riu)du+ I -rrdu. (C19) 

It also shares the properties 1, 2, and 3. 

Observe that the left-hand side of (C.16) is an odd function in the annulus 0° \ D^^g. Hence 
the right-hand side is odd as well. Applying this to z = id2 we obtain the equation 

/ (f^o{u) + N-^^^) du + N-%{id2) + N-''/^Ao{id2) + c, 

[ (fxoiu) + N-^^^) du + N-'fo{-id2) + N-'/^Ao{-id2) + ci 
A1U72 V 2fio{u)/ 



IJ 71 U72 

The functions iJ,o{z) and U-^{z) are even, hence 



+ 0{N-^). 
((7.20) 



mn TTi 



Ho(u) + N-^ '-^]du= Uo(u)+iV-i -A]du , (C.21) 

^ 2mo(u); 4 V ^ ' 2/xo(«)y N 2N' ^ ' 

with the term — ^ — coming from regularization. Indeed, due to (C.13), as z — 00, Imz < 0, 

l^o{z) = -fisz^ - fiiz - fi_iz~'^ +Jio{z), Jlo{z) = 0{z~^), (C.22) 
and we define the regularization of ^ IJ'o{u)du as 

/ iJ,o{u) du = —a^ + —a'^ + H-ilna + / Jio(u) du + / iio{u)du, (C.23) 
J 13 42 7„ 

where z = —id2 and a is any point on 73. The integral is taken over the part of 73 from a to 
00 (in the lower half-plane). Similarly, as z — > 00, Imz < 0, 

= 1-1 + r{z), r{z) = 0{z-% (C.24) 
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and we define 

/ V4 du = + / r(u) du + / V4 du . (C.25) 

4 2/"o(«) 2 J, ^ > A 2/io(n) ^ ^ 

We take the main branch of logarithm for Ina in (C.23,25). The remainder functions ]jio{z) and 
r{z) are even. Since 

ln(zc?2) = ln(-zd2) + m (C.26) 

and ii-i = - jj, 

L +^"^) = T^^^^)' + f^^^^)^ + ('^-^ - 2^) 



W~2N' 



hence (C.21) is proved. 

We obtain from (C.20,21) that 



2i2o{u) J 2N AN 2 



^_5/3 Ao(id2^AoH^ + 0(7V-2). 



(C.27) 



For ix G 7, 

-I- A/"~i 

2//o(?i) 

hence 



AioH+A^-'T^^^= [-cZ(n) + Ar-iC/i(t^)]'/' + 0(Ar-2n-2), (C.28) 



[ (^^^{u) + N-'^^^)du= [ [-d{u)+N-^U\u)Y^^du + OiN-^). (C.29) 
J J V 2^o[u)J 

Deforming the contour of integration we obtain that 

r /'oo poo 

[-d{u)+N-^U^{u)Y'^ du = -2 [-d{u)+N-^U^{u)] ^ du = -2 ij,%u)du. (C.30) 

Thus, (C.27) reduces to 

^ H\u)du + ^ + + N-^a + N-^'H + 0{N-% (C.31) 

where 

_ fo{id2) + fo{-id2) Ao{id2) + Ao{-id2) (nr,n\ 

a — 2 ' " — 2 ■ (.o.^^j 

Let us evaluate a and 6. 
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Prom (C.15), 



2 a?2(-^) 2 a02(2;) 2 fio{z) - a1^{z) 



(C.33) 



Prom (4.4) and (C.4), 



ai2(^) ai2(^) a"i2{z) 



{CM) 



hence 

Thus, from (C.33) (and the inequahty < Imfo(2;) < f ), 

2 ^' a?2(^) 



fo{z) + fo(-2) = ^ - iV-2/3_J_ + O(iV-i), (C.36) 



and 

Let us cvahiate Ao(-z) + Ao(— -z). 
Prom (C.6'), 

Ao(^) = Ao(^) + Ao(z), (C.38) 



where 



The function Ao(-2;) is odd, hence Ao{z) + Ao{-z) = 0. Por Ao{z) we have that modulo 0{N ^/^)- 
terms, 

Ao(2;) + Ao(-z) = ( -Q — + — 







= -2c0 



-2c0 



-2c° 



-/xg + 2zMoa0i + (aOi)2 + (a?i)2 

i/xp + g^i 

-a^saOi + 2i/xoaOi + {^^f 



(C.40) 



2a^i (i/xo + ^21 ) / a^i (^) ^12 (^) 



Thus, 



= -iV"^/^ Q +0{N-^), z = id2. (C.41) 



Combining (C.37) and (C.41) we obtain that 



N-^a + N-^/^b = -— + 0{N-^), ((7.42) 
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hence by (C.31), 

c,= I ^'^iu)dn + ^+0{N-% 

so that by (C.15), 

1 Trin 

^ H%u)du - _ Ini^o + — + 0(iV-2). (C.43) 

-0 

From (C.12) and (C.43) we obtain (5.14). Proposition 5.2 is proved. It remains to prove (C.14). 

Proof of (C.14)- To simphfy notations we will drop super- and subscripts '0'. In the proof we 
will neglect 0{N~^) terms. We will prove first that 

t' = f' + f (C.44) 

and then we will check (C.14) at infinity. Prom (C.15), 

, - a'n a'i2 , «ii«i2 - aiiai2 ai2(/x/x' - ana'^i) + a'i2(an - /x^) 

2(r + r ) = 1 = r . 

/x-oii ai2 /xai2 ai2yw(yw - On) 

Since 

/x^ = On + 012021, 2/x/i' = 2ana[i + o'i202i + 012O21, 

we obtain that 

4(f' -|- f ) — ^12(012021 + 012021) — 2oi2ai202i _ 012O21 — 012021 

oi2/x(/i - On) /x(/x-aii) 

which coincides with 4r', see (5.5). This proves (C.44). As z — cxd, Imz > 0, 

r{z) = 0{z-^y, f{z) = Unz-UnRl + 0{z-^)- f'{z) = -\z-^ + r{z), f{z) = 0{z-^), 

hence 

1 r ^ 

f (z) = — - In^; + / r{u)du. 

2 J 00 

This implies that 



lim t{z) = lim 

z— >oo z— >oo 

Equation (C.14) is proved. 



f{z)+f{z) + hnRl 



0, Im^; > 0. 
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Appendix D. Proof of Lemma 6.2 

Proof of (6.27). Let us rewrite (6.24) as 



{D.l) 



where 



1 1 
-1 1 



] = ( Ql2 
/ V /^^(^) J 
1 

_Ml£)+a2j>) 



1 

aii(^) 







1 1 
-1 1 



a;,(z) J 



{D.2) 



and 



1/2 



Prom (1.68), 



Ui^z) = {U%z)Y" = {-d{z) + N-'U\z)y'' = n{z) + N-'^ ^ + o{N-'\z\-% 



il U\z) 



2 i^iz) 



U\z) = {a^,J{z)-aUz)^^^'^"^ 



{D.4) 



therefore we can replace iJ,'^{z) for ij,(z) in formula (D.2) for T'^(z), with an error term of the order 
of N~^\z\~^ which we will omit. This gives 



T%z) 



3° (z)\^^^ /I 1 

,0 ^.^\ 1/2 / 1 «?2(^) 



H{z) - al^jz ) 
H{z) 

To{z)e-^^'^''' 



1 



where 

hence by (D.l), 
Prom (D.3), 

where 

Prom (D.4), 



f(z) = - In 
^ ' 2 



M(z)-aii(^) 



/i(z) -a?i(z) 



(D.5) 



a?2(z) 



{D.6) 



J oo 

= ii\u)du. 



{D.l) 
{D.8) 



/ iJ,^{u)du= ii{u)du + N-^f{z) + 0{N-^\z\-^) =^{z) + N-^f{z) + 0{N-^\z\-^), {D.9) 

J oo J oo 
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where 

^z) = 1 r du. {D.IO) 

Thus, modulo terms of the order of A'^~^|z|~^, 

*wkb(^) = Co^o(z)e-[^«(^'+^(^)+^(^)+^>^ {D.U) 
Prom (5.5), (D.6,10) we obtain that 

f'{z)+T'{z) = T'{z) {D.12) 
(cf. Appendix C above). In addition, as z — > oo, 

f{z) = ilnz- i Initio +0(z-2); f'{z) = -\z-^- r' (z) = 0{z-'). 

Hence, integrating (D.12) we obtain that 

f{z)+fiz)=T{z)-^lnRl {D.13) 

where integral (D.IO) is regularized as 

1 r fuHu) 1 



Thus, we obtain that, modulo terms of the order of AT \z 



u 

lui-i 



du. 



*wkb(^) = CoTo(z)e-[^«(^)+-(^)+^^l'^3 = ^^kb(^), C^=C^-\ ln<, 
which was stated. 

Proof of (6.28). Let z G Tq. For the sake of definiteness, assume that Imz > 0, so that 

*TP {z) = C^W{z)Y^{w{z)), 

where 



From Proposition 6.3, 

= Ai {N'^/^z) 



— TT + £ < arg Z < TT — £, 



and 



27r 

TT + £ < arg z ^ < TT — £. 

o 
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Thus, 



_ 1 /^-l/4 \ / e-2JV^^/V3 ^2Nz^'^/Z\ 



TT 

— + £ < arg z < TT — £. 



Prom (6.12), 



11 (ujdu = z, \z), 



and from (6.10), 



hence 



(.) = C,W{z) (1 + 0(iV-)) ^ - l^^.,^, j 



TT 

— — + £ < arg t«(z) < TT — £. 
o 



Observe that if z G r° and Im z > then 



TT 

— — + £ < arg w{z) < TT — £, 
o 



hence (D.20) holds. Prom (1.79), 



J V a?,(z) a?,(z) 



hence from (D.20), 



*TP (z) = (1 + O(iV-i)) ^r^(^)£;(iV^^(^)) = (1 + O(iV-i)) *wkb(^), 



if Co = Lemma 6.2 is proven. 
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Appendix E. Proof of Lemma 4.2 

Proof. Let i? > be a big number, independent on N, which will be chosen later. Consider 
two cases: (1) \z\ < RN'^/^, and (2) \z\ > RN'^/^. In the case (1) detW{z) / by (4.63). 
Consider the case (2). First we notice that by (4.48), 



z G (-J^i)UO°UJ^i. 



Therefore by (4.67), 

011(2;) 
021(2;) 

By (4.4), 



-a22{z 

v2 



0{N- 



-1/3, 



012(2) 
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,2_^2 + 0(iV-2/3), 



aUz) = -aUz) = + 0{N-'/'z), a%{z) = [Rlf/^gz^ + 0{N-''% 

al,{z) = -{R^S/^gz^ + 0{N-^'^). 
By (4.66) this gives the matrix elements Wij{z) of W{z) as 



■wii{z) 



i-z^ + {Riy/'gz' + 0{N-'/% 



wMz) = W2i{z) = ^ + 0{N-'r^z), 



VJ22{z) 



gz 



,2 2 + (^0)1/2^^2^0^^-2/3)^ 



and these estimates are valied for all z G (— i^i) U U J7i. Observe that by (1.45,53), 



Also, 

for all Izl > i?Ar-i/3. Hence 



= 0{R-^), 



W{z) = 



gz 



y^Z^- z'^ +Zo 



y/zQ - Z^ + Zo 



+ 0{R-^) 



for all z e i-ni) un^UQi and \z\ > RN-^/^. Since 
^/zq- z'^ + zo z 



det 



^/zq- z^ + Zo 



2a/z, 



2^ + 20 /0, 



(E.l) 



{E.2) 



(E.S) 



{E.3) 



{EA) 



{E.5) 



for z e (-l^i) U f]" U J7i, we obtain that det W{z) 7^ for sufficiently large RiizE (-O1) U 0° U fii 
and l^l > RN~^^^. Lemma 4.2 is proved. 
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Appendix F. Alternative Forms for *^(^;) 

Recall that each of the functions ^{z) and Co{^) defines an analytic change of variable in the 
rectangle = U $7° U f^i, and according to (4.84) we have that 

C(z) = CoW + 0(Ar-4/3), zen,. (F.l) 

Let us use the function Co (-2), instead of the function the right hand side of (5.10). This 

will lead us to the matrix function 

' CV{z)^'' (n^/^Co{z)) , lmz>0, 

*CP {z) = { ) ( (^-2) 

\^CV{z)^'^ (^N^/^Co{z)j , lmz<0. 

It is clear that the functions '^^p (z) and (z) have exactly the same Stokes matrices. Taking 
also into account that, by (F.l), 

g-i(|iVf'^(z)+ivi/3yC(z))a3gi(|iVC|(^)+ivi/32/Co(2))<^3 = 7 + 0(Ar-V3)^ 

we arrive to the asymptotic formula, 

^Cp{zm%p{z)]-'=I + 0{N-'/'), zeOP. (F.3) 

The last equation together with (8.21) yield the following representation for the solution ^n{z) in 
the domain 0°, 

^r^{z) = (l + 0(iV- 1/3)) $0 p (^), ^ e a^, ^FA) 
Of course, equation (F.4) has a bigger error term then a similar equation 

^'^(z) = (/ + 0(iV-i))*cp(^), ze^\ (F.5) 

which is given by Theorem 1.2. Sometimes, however, a certain advantage can be achieved by 
using (F.4). Indeed, the function C,o{z), unlike the function C,{z), is given by an elementary explicit 
formula (4.49). It is also worth noticing that estimate (F.4) can be improved as follows. 
Consider the matrix ratio, 

R{z)^^cv{z)[^%P {z)]-\ 

The function R{z) does not have a jump on the interval [— and hence it is analytic in ^2°. 
Therefore, 

^(,)_/- -l>cp(»)fe(»)l ^^„„, 
Jqqo u- z 2m 

At the same time, on the boundary of the domain the function ^'cp (z) can be replaced by 
*wkb(-2) with an error of the order of N~^, 

*cp {z) = {I + 0{N-^)) *wkb(2), z G (F.6) 

and therefore, 

^ /■ (/ + 0(Ar-i))^wKB(n)[^^p(n)]-i du 
Jqqo u- z 2m' 
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In addition, it follows from (F.6) and (F.3) that 

^WKB{n)[^% (u)]-' =1 + 0(iV-i/3), z G dn°, (F.8) 
and hence equation (F.7) can be rewritten as 

/■ ^wKB(n)[^°cp(")]"^ I ^ ^ (F^) 

^'~Jano u-z 27rz V^dist(2;5J^o)y'- ^ ' 

The last equation, together with the standard arguments based on the flexibility in the choice of 
the basic geometric parameters di and 6,2, yields the following improvement of estimate (F.4) 

*„(^) = (7 + 0(Ar-i))*i,p(z), (F.lO) 

where 

JdQ.0 u-z 2m J 

The only non-explicit element, apart of course from the Painlevc functions, which is left in formulae 
(F.IO) and (F.ll) is parameter a in definition (4.82) of the variable y. To determine it one needs 
equation of periods (4.78). However, the terms of order N~^/^ can be extracted from (F.lO)-(F.ll) 
in the very explicit form. Indeed, by a strightforward though a little bit involved calculations one 
can see that 

*wkb(^)[*°cp i^r' = /-iV-V3^^M= (zas - HR^J'/'a^) +0{N-y% z G dQ^, (F.12) 

V ^0 ^ 

where 



2Coo(^) gj, \ V?-4) U^^ju^-Zl Coo{z) 

By a direct calculation (this time much less involved) we check that ri{z) has no pole at z = 0, i.e. 
it is an analytic function in O*'. Therefore, equation (F.ll) yields the estimate, 

*J,P (.) =(^- N-'''^= H - 2^{RlY'^a,) + 0{N-^'^)^ (.). (F.13) 

Further improvement of (F.13) could be done. However, when we go beyond the order N~^/^ 
we have to take into account the parameter a and hence representation (F.IO) will not give us any 
advantages comparing with our basic formula (8.21). 

In principle, it is possible to avoid the use of period equation (4.78). To this end, let us 
define the approximate solution '^^{z) of the Riemann-Hilbert problem by relation (8.1) where 
^CP {z) is replaced by ^'^p {z) and a in the definition of y is put to zero. Then for the ratio 
Xn{z) = '^niz)[^niz)]~^ wiU arrive to the Riemann-Hilbert problem whose jump matrix G{z) = 
[Xn-{z)]~^Xn+{z) is given explicitely in terms of *°(2;) and satisfies the uniform estimates (cf. 
(8.11) and (8.12)) 

ll^-G'(z)|U,(^)nL^(^)=0(Ar-V3). 
This would mean the replacement of the final estimate (8.21) by the estimate 



*„(.) = {i + o ij^TjsjYTW) ' ' *° ^^'^^^ 



70 



Theoretically, this estimate can be improved by iterating the corresponding singular integral equa- 
tion (cf. [DZ3]). In particular, the first iteration leads to the function ^'^p (z) satisfying (F.13). 
The second iteration, however, looks extremely cumbersome. In fact, our basic construction involv- 
ing the nontrivial change of variable function ({z) can be thought of as a way to bring the second 
iteration of the singular integral operator with the weight / — G{z) to a compact form. 

Let us finally discuss the dependence of Co(-2^) on the parameter y, Co{z) = Co{z',y)- To that 
end put 

M,{z, Co{z; y), y) ^ {z, Co{z; y), y)e(^f ^o(^^^)+^'''^^«o(^^^)-^t)'^3 ^ 
z E {e < arg z < ir — e}U {— tt + e < arg z < — e}, 



(F.15) 



and 

M2{z, Co{z; y),y) = {z, Co{z; y), y)5ije(^f ^o(^^^)+^'''^^^°(^^^')-^^)-3 , 



(F.16) 



r TT TT , T ^1 

z E \ —— + e < arg z < — e, ±Im z > > 
[the matrices S^ l are the Stokes matrices from (7.3)]. It is not difficult to see that 

Mj{z,(:o{z;y),y) = (l + 0{N-"^)) M,{z,Co{z;Q),y). {F.17) 
This implies that we can redefine {z) as 



2; G {e < arg 2; < tt — e} U {— tt + e < arg z < — e}, 

and 

r TT TT , T ^1 

z G | — — -I- e < argz < — — e, ±lmz>Oj, 
and still have (F.4). Observe that by (4.49,44), 



(F.18) 



(F.19) 



Co{z;0)=\^l£^^z^,-u^du 



1/3 



(F.20) 



Formulae (F.18)-(F.20) are of the type that would appear if for the analysis of the Riemann-Hilbert 

problem (1.35-38) we use the nonlinear steepest descent approach similar to the method used in 
[BDJ] for the orthogonal polynomials on the circle. The advantage of such scheme would be the 
apearence of the Painleve Riemann-Hilbert problem in a completely deductive way, i.e. without use 
of any prior information about the structure of the asymptotics of the coefficients i?° . The price, 
however, would be the much more subtle structure of the relevant contour deformation and much 
more complicated analysis of the all error terms than in our approach. In our approach, we take 
the full advantage of the prior heuristic study of the Freud equation (1.24) which leads us to the 
choice of the approximate solution, whose justification only needs a very "light" Riemann-Hilbert 
analysis. Moreover, we obtain the error term up to the order 0{N~^). To get this accuracy in 
the framework of the [BDJ] scheme would in fact also require, a posterior, a WKB-type analysis 
of the associated Lax pair (cf. [DZ3]). In our method, we use the WKB formulae from the very 
beginning. 
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